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[Text]  Abstract 

A  parallel  matrix  multiplier  based  on  only  a  few  imaging  systems  is  presented. 

Many  computation  problems  in  science  and  engineering  eventually  become 
algebraic  operations.  Hence,  the  development  of  special  devices  to  solve 
complicated  algebraic  operations  has  always  been  desirable. 

This  paper  presents  a  parallel  matrix  multiplier  which  involves  only  a  few 
imaging  systems.  It  is  based  on  the  corresponding  decomposition  of  optical 
parallel  vector-vector  out-product  multiplier  and  matrix  multiplier.  It  can  be 
used  to  solve  real  analog  multiplication  of  matrices,  analog  two-dimensional 
transforms,  digital  precision  matrix  multiplication,  and  digital  precision 
two-dimensional  transforms.  It  takes  full  advantage  of  the  parallel  nature  of 
optics  in  all  operations.  The  multiplication  of  three  matrices  or  a  two- 
dimensional  transform  takes  place  within  one  clock  cycle  of  the  system. 

I.  Structure  of  the  Parallel  Matrix  Multiplier 

1.1  Multiplying  Two  Matrices 

Matrix  multiplication  can  be  decomposed  into  a  number  of  vector  inner-products 
and  outer-products  (or  tensor  products). 

Figure  1  shows  a  novel  optical  parallel  matrix  multiplier  based  on  vector- 
vector  out-products .  The  system  employs  n  point  light  sources  on  the  focal 
plane  of  its  spherical  lens  system.  In  the  figure,  SL  and  CL  are  spherical  and 
cylindrical  lenses,  respectively.  Their  focal  lengths  obey  the  following 
equations : 
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fsLx  ^  ^CLx  ^  f  CL2  ~  fcL^  ^  f  * 
fa,  =  3//4 

In  the  figure,  A  and  B  are 
spatial  light  modulators 
(SLMs)  used  to  input  the 
matrices .  They  are  called 
input  encoding  plates .  The 
transmittance  of  each  element 
is  proportional  to  the  value 
of  the  corresponding  matrix 
element.  Therefore,  we  have  to 
assume  that  every  matrix  element  is  a  positive  real  number.  As  for  different 
situations  involving  numbers  that  are  not  positive  real  numbers,  there  are  a 
variety  of  methods  available.  Even  a  complex  matrix  can  be  coded  to  become  a 
real  matrix  before  input.  In  the  system  shown  in  Figure  1,  addition  and 
multiplication  are  done  in  the  same  manner  as  with  an  inner-product  matrix 
multiplier.  Figure  1  shows  the  principle  of  the  optics.  Obviously,  it  is  an 
orthogonal  imaging  nonspherical  optical  system. 

If  the  light  sources  are  ideal  point  sources  and  all  lenses  are  ideal  lenses 
(without  taking  aberration  into  account  and  cylindrical  lenses  only  form 
images  in  one  direction)  ,  then  matrix  code  plates  A  and  B  are  placed  as  shown 
in  Figure  1.  If  the  point  source  intensities  are  Io5(X-Xo)5(Y-Yo)  ,  Ii5(X-Xi) 

5(Y-Yi),  l2^ (X-X2)5(Y-Y2)  ,  . (X-Xn-i)5(Y-Yn-i)  ,  and  the  transmittances 

of  encoding  plate  A  and  B  are  (A,-j)  and  (Bij),  then  it  is  apparent  that  there 
is  the  following  intensity  distribution  on  plane  P: 

c'li  =  gY,  AiAjIj, 

k-0 

If  1(5  -=  Ii  -  .  “  In_i  -  I,  we  have: 

Gij  =  C'ij/Gl  =  .  (1) 

This  is  the  product  of  matrix  A  and  matrix  B.  Of  course,  transmittance  is 
always  less  than  1.  In  reality,  it  is  necessary  to  multiply  by  a  proportion¬ 
ality  factor.  As  for  the  rigorous  derivation  of  equation  (1),  it  can  be 
obtained  by  using  a  transfer  function  technique. 

1.2  Multiplication  of  Three  Matrices 

If  we  need  to  multiply  three  matrices,  then  the  system  shown  in  Figure  2  may 
be  used,  the  difference  is  the  light  source.  In  Figure  2,  n  x  n  point  sources 
(assuming  matrix  B  is  n  x  n)  are  used,  if  the  three  matrices  to  be  multiplied 
are  (A^j),  (Bij)  1  and  (C^j),  then  A  and  C  are  placed  in  the  manner  shown  in 
Figure  2.  The  light  intensity  of  each  light  source  is  proportional  to  the 
matrix  B,-j  (i,j  “0,  1,  . ,  n-1) ,  and  in  their  relative  positions  are  also 
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Figure  1,  Optical  Matrix  Multiplier  Based  on 
N  Vector-Vector  Outer  Products 
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shown  in  Figure  2.  It  is  apparent  that  the  light  intensity  distribution  on 
plane  P  is: 


(2) 
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where  the  s3naibol  repre¬ 

sents  "proportional  to . "  This 
is  the  product  of  the  three 
matrices.  Please  note  that  the 
order  of  the  three  matrices  in 
equation  (2)  is  different  from 
that  in  Figure  2 .  This 
arrangement  provides  some 
convenience  for  the 
multiplication  of  more  matrices 


Figure 
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Optical  Three-Matrix  Multiplier 
Based  on  Vector-Vector  Outer  Product 


II.  All-Parallel  Optical  Multiplier  Technique  and  Its  Optical  Two-Dimensional 
Transform 


Because  light  intensity  is  always  positive,  all  matrix  elements  described  in 
the  previous  section  are  positive.  In  the  multiplication  of  matrices,  some 
elements  may  be  real  or  complex  nvimbers.  This  requires  us  to  take  certain 
measures . 

This  paper  presents  a  technique  which  involves  a  combination  of  polarization 
encoding  and  light  intensity  encoding  to  realize  optical  two-dimensional 
transform. 

2.1  Expression  and  Operation  of  Intensity/Folarizatlon  Hybrid  Encoding  for 
Real  Numbers 


(1)  Expressing  a  real  number 


To  express  a  real  number  by  intensity  and  polarization,  the  magnitude  of  light 
intensity  can  be  used  to  express  the  absolute  value  and  two  different 
directions  of  polarization  can  be  used  to  represent  different  signs.  Based  on 
this  concept,  an  arbitrary  real  number  may  be  expressed  as  follows: 

M  -  |M|[1  +  sign  (M)]/2(+l)  +  |M|[1  -  sign(M) ]/2(-l) 


where 


8ign(H) 


+1  «>o 
0  Af=0 
-1  M<0 


(3) 


Optically,  ex  and  ey  are  two  independent  orthogonal  polarization  directions. 
It  is  further  assumed  that  ex  corresponds  to  (+1)  and  ey  to  (-1).  Then,  the 
light  beam  representing  the  number  M  has  the  following  two  polarization 
components  {Ix,  ly) : 
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(4) 


/j,  =  |/^|  [i+sifim(«)]/2 

]ly  =  \M\  {l-sigtn(W)  ]  /2 

TllUS  f  M  “  lyCy 

and 

Ixly  =  0  (5) 

Obviously,  this  expression  is  unique.  In  order  to  meet  the  need  in  practice, 
the  two  independent  polarization  components  of  the  light  beam  expressing  M  can 
be  expressed  in  a  more  general  manner  as : 

fj,  =  |W|  [l  +  sism(A0]/2  +  C  /g) 

|Jy=  |W|[l-sign(Af)]/2+C 
C  >  0 

In  this  case,  the  expression  is  not  unique. 

(2)  Arlthinietlc  operation 

The  multiplication  of  two  real  numbers  can  be  expressed  by  the  method  shown  in 
equations  (4) and  (5),  i.e. 


M  “  ImxSx  +  lMy®y 
N  "  Inx®x  ^My®y 

J^=  |W|  [1 + sign  (Af)]/2 
=  jw|[l-sigcn(«)]/2 

=  |JV|[l+sifim(W)3/2 
|i/|  [1-sign  (iff)  ]/2 

MN  -  |M||N|tl  +  sign(M)  sign(N)]/2ex 
+|M||N|[1  -  sign(M)  sign(N)]/2ey 
-  IMNjil  +  sign(MN)]/2ex  +  |MN| [1  -  sign(MN) ]/2ey  (7) 

The  multiplication  of  two  real  numbers  optically  expressed  by  equations  (4) 
and  (5)  can  be  realized  by  shining  two  polarized  light  beams  on  two 
transmission  plates.  If  the  transmittances  are  rji,  tf2,  then 

^  “  ’?1’72[^0x®X  3ioy®y]  (^ox^oy  “  0) 

In  order  to  express  the  sign  of  operation  of  the  two  nximbers,  it  is  possible 
to  stick  a  layer  of  polarizer  on  the  transmission  plate.  Transmittance  and 
polarization  angle  together  can  express  a  number  M  with  a  sign.  When  M  <  0, 
then  a  layer  of  polarizer  can  be  used  to  rotate  the  direction  of  polarization 
of  the  incident  light  by  90®.  When  no  polarizer  is  required  and  the 

direction  of  polarization  remains  unchanged.  In  the  meantime,  iji  -  |M| .  The 
intensity  of  the  incident  beam  represents  the  value  of  N  and  its  direction  of 
polarization  represents  the  sign  of  N.  Then,  the  output  is  (M*N),  which  in 
equation  (7)  is  the  product  of  two  real  numbers. 
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If  addition  or  subtraction  is  required,  then  it  is  necessary  to  use  the  method 
shown  in  equation  (6).  In  this  case,  the  expression  is  no  longer  unique.  For 
example,  the  subtraction  of  two  greater  than  0  real  numbers  A  and  B,  i.e., 

C  “  A  ■“  B,  can  be  expressed  as: 

C  *  Iax®x  “*■  ^By^y 

In  this  case,  subtraction  becomes  adding  a  positive  value  to  the  polarization 
component  which  represents  the  negative  value. 

Obviously,  the  number  must  be  transformed  after  addition  or  subtraction  before 
multiplication  can  be  carried  out  again.  Nevertheless,  it  is  very  difficult  to 
do  so  directly  by  optical  means.  Hence,  it  is  usually  done  electrically. 

2.2  Two-Dimensional  Transform  of  Intensity/Polarization  Hybrid  Code 

Any  two-dimensional  linear  transform  in  the  complex  domain  can  be  reduced  to 
four  independent  real  number  transforms.  Therefore,  this  section  only  deals 
with  two-dimensional  linear  transforms  in  the  real  number  domain.  (This 
remains  true  unless  specifically  expressed  otherwise.) 

Any  two-dimensional  transform  can  be  expressed  as  follows: 

E  -  UFV 

where  F  is  the  matrix  to  be  transformed,  U  and  V  are  the  transformation  core, 
and  E  is  the  end  result.  If  F^j  is  expressed  by  the  method  described  earlier. 
|U|-j|  and  |V^j|  are  expressed  in  terms  of  encoding  plate  transmittance  and  sign 
(U,*j)  and  sign  (Vf j)  are  expressed  by  the  polarizer  sticker  method.  The  optical 
system  to  realize  this  transformation  is  shown  in  Figure  2.  Its  light  source 
is: 


Fij  ■=  |F,j|fi(x-Xi)5(y-yi){  [1  +  sign(Fjj)  ]/2ex  +  [1  -  sign(F,j)  ]/2ey) 

The  encoding  plates  are: 

Uij  =  |Uij|fi(x-Xi)5(y-yi){  [1  +  sign(Ujj)  ]/2ex  +  [1  -  sign(U,-j)  ]/2ey) 

^ij  “  iv,jjS(x-x,-)5(y-yi){  [1  +  sign(Vij) ]/2ex  +  [1  -  sign(V,j) ]/2ey) 

For  a  specific  type  of  two-dimensional  transform,  Ujj  and  Vf j  remain  unchanged. 
Hence,  they  can  be  constructed  as  fixed  encoding  plates.  Thus,  the  system  does 
not  require  an  SIM  as  the  input;  it  only  needs  two  light  arrays  that  can  be 
modulated.  If  F|j  is  always  greater  than  0,  it  is  even  simpler.  It  only  needs 
one  light  source  array,  which  is  extremely  convenient. 

In  Figure  1,  all  we  have  to  do  is  to  subtract  the  two  polarization  components 
received  correspondingly.  This  only  requires  equation  (6),  and  the  result  is 
apparent.  The  subtraction  of  two  polarization  components  can  be  completed  on 
an  electronic  circuit .  There  is  no  need  to  do  it  optically. 
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In  conclusion,  the  principle  of  an  optical  parallel  vector-vector  outer- 
product  multiplier  is  presented.  It  is  very  simple  and  can  easily  be  expanded 
to  larger  dimensions.  By  using  common  optical  components,  it  is  possible  to 
obtain  the  outer-product  of  10*  x  10*  vectors  (in  analog).  The  system  can  be 
used  to  realize  multiplication  involving  two  and  three  (positive  real  number) 
matrices  (in  analog).  Experimental  results  to  verify  this  theory  will  be 
published  separately.  (See  next  issue.) 

On  the  basis  of  hardware  research  on  an  analog  matrix  multiplier,  the 
algorithm  is  discussed  in  this  paper.  It  expands  the  range  of  applications  to 
the  multiplication  of  real  number  matrices.  Furthermore,  a  hybrid  spatial 
intensity/polarization  encoding  method  is  used  to  realize  the  multiplication 
of  two  to  three  digital  (mixed-base)  matrices  to  improve  the  accuracy  of 
operation.  Compared  to  the  two  existing  all-parallel  matrix  multipliers  (based 
on  inner-product  and  Fourier  transform) ,  this  system  has  a  distinct  advantage 
in  performing  two-dimensional  linear  transform  in  the  real  number  domain.  It 
only  needs  a  light  array  to  represent  the  pattern.  The  transformation  core  is 
a  number  of  fixed  encoding  plates.  The  result  is  obtained  in  a  single  step.  It 
does  not  require  an  SLM.  (This  is  because  matrix  B,  in  the  multiplication  of 
matrices  A,  B,  and  C,  is  expressed  by  a  light  array  in  this  approach.  In  other 
approaches,  matrix  A  is  expressed  by  a  llg^t  source  Instead.)  Therefore,  use 
of  this  triple  matrix  multiplier  system  should  be  promoted  for  two-dimensional 
transforms . 


Experimental  Research  on  Optical  Parallel  Multiple  Matrix  Multiplier 
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[Text]  Abstract 

The  experimental  method  and  results  of  a  parallel  optical  matrix  multiplier 
are  described.  The  system  takes  full  advantage  of  optical  parallelism  and  is 
capable  of  handling  the  multiplication  of  three  matrices  of  a  two-dimensional 
transform  within  one  clock  cycle.  The  sources  of  error  and  ways  to  eliminate 
them  are  discussed.  Results  of  multiplications  involving  two  and  three 
matrices  (analog  operation)  are  given.  The  mean  error  of  two-matrix  multipli¬ 
cation  is  1.2  percent  and  that  of  three-matrix  multiplication  is  1.23  percent. 

Currently,  matrix  operations  are  done  serially  in  computers.  Hence,  the 
inherent  parallelism  of  matrices  is  totally  lost.^'^  Optical  computing  is 
extremely  parallel  in  nature.  Therefore,  researchers  are  actively  seeking 
optical  means  to  realize  parallel  matrix  operation.  To  date,  there  are  many 
optical  matrix  multipliers  which  have  been  briefly  discussed  in  reference  [3] . 
A  parallel  matrix  multiplier  involving  only  a  few  imaging  elements  was 
introduced  in  reference  [3].  Multiplication  of  real  number  matrices  is 
realized  based  on  optical  parallel  vector-vector  outer  product  and  the 
corresponding  decomposition  of  matrix  multiplication.  This  paper  describes  the 
experimental  process  and  results  of  this  parallel  matrix  multiplier.  It  makes 
full  use  of  optical  parallelism.  The  multiplication  of  three  matrices  or  a 
two-dimensional  transform  takes  only  one  clock  cycle. 

Principle  of  the  All-Parallel  Matrix  Multiplier 

Figure  1  shows  the  vector-vector  outer  product  based  optical  parallel  matrix 
multiplier  presented  in  reference  [3].  The  system  employs  n  x  n  point  sources 
placed  on  the  focal  plane  of  a  spherical  lens  system  for  Illumination. 
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In  addition,  it  represents  the  matrix  B.  In  the  figure,  SL  and  CL  are 
spherical  and  cylindrical  lenses,  respectively.  Their  focal  lengths  follow  the 
equations  below: 

fsi^  -  fell  ~  ^  f 'fexj  “  3//^ 


The  SLMs  used  to 
input  the  other 
two  matrices  are 
A  and  C.  In  our 
preliminary 
experiment ,  they 
were  replaced  by 
input  coded 
plates.  (In  a 
two-dimensional 
transform,  these 
two  coded  plates 
represent  the 
core  of  trans¬ 
formation  . )  The 
optical  trans¬ 
mittance  of  each 
element  is  pro¬ 
portional  to  the 
value  of  the 
corresponding 
matrix  element. 

To  this  end,  we 
have  to  assume 
that  all  matrix 
elements  are 

positive  real  numbers.  When  they  are  not  positive  real  numbers,  there  are  a 
number  of  schemes  available.  Even  when  they  are  complex  numbers,  a  matrix  can 
be  coded  to  become  a  positive  real  number  matrix*  prior  to  processing.  In  the 
system  shown  in  Figure  1,  optical  addition  and  multiplication  are  accomplished 
in  the  same  manner  as  calculating  the  inner  product  of  matrices.  Figure  1(a) 
Illustrates  the  principle  of  the  optical  system  and  Figure  1(b)  shows  the 
paths  of  light  rays.  Apparently,  it  is  a  nonspherical  orthogonal  imaging 
optical  system. 


Figure  1.  (a)  Optical  Scheme  of  the  Three-Matrix 
Multiplication  System 
(b)  Paths  of  the  Light  Rays 


If  the  light  sources  are  ideal  point  sources  and  all  lenses  are  ideal  imaging 
lenses  (without  considering  any  aberrations,  and  the  cylindrical  lens  forms  an 
image  only  in  one  direction) ,  then  coded  plates  A  and  C  are  placed  according 
to  the  way  shown  in  Figure  1(a).  If  the  li^t  intensities  of  the  point  sources 
are  Bip,,  the  transmittance  values  for  coded  plates  A  and  C  are  (A{j)  and  (Cfj), 
then  the  receiver  plane  should  have  the  following  intensity  distribution: 


=  D^G  = 


(1) 
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This  is  the  multiplication  of  matrices  A,  B,  and  C.  Of  course,  transmittance 
is  always  less  than  1.  In  reality,  it  is  always  necessary  to  multiply  by  a 
proportionality  constant. 

For  the  multiplication  of  two  matrices,  all  we  have  to  do  is  to  change  matrix 
B  to  a  unity  matrix,  i.e.,  to  change  the  light  source  to  n  point  sources. 

Analysis  of  the  Experimental  All-Parallel  Optical  Matrix  Multiplier  System 

The  major  sources  or  error  for  the  matrix  multiplier  described  above  include: 
(1)  error  due  to  diffraction  of  small  elements  on  the  coded  plates,  (2)  error 
due  to  nonunifora  optical  angular  distribution  of  different  point  sources, 
(3)  error  from  background  light,  (4)  error  caused  by  distortion  of  coded 
plates,  and  (5)  error  due  to  deviation  of  the  actual  optical  system  from  an 
ideal  system.  When  the  size  of  the  elements  of  coded  plates  are  large,  (2), 
(3),  and  (4)  have  larger  impact  on  the  error.  It  is  easy  to  reduce  the  last 
one  to  meet  the  accuracy  requirement. 

1.  Relation  between  size  of  optical  elements  and  coded  plates 

Let  us  assume  that  the  point  light  sources  are  Sx  and  Sy  from  center  to  center 
and  the  elements  on  plates  A  and  B  are  d^x,  d^y,  d2x  and  d2y  from  center  to 
center,  respectively,  and  each  element  of  receiver  is  a  small  block  hx  x  hy  in 
size.  If  the  system  is  an  ideal  optical  system,  then  based  on  optical  wave 
theory,  the  resolution  light  is 


hxi  «  1.22A/NAj 

hyi  «  1.22A/NA, 

Let 

Then 

hx  “  hy  —  h>2h: 

It  is  obvious  that 

Sx  “  Sy  “  dix  “d2x  “  diy  “  ^2y  “  ^ 

The  aperture  of  the  system  should  be  determined  by  the  outer  frame  of  A  and  B. 
Hence,  there  are  specific  requirements  on  minimum  aperture  for  every  lens. 

In  order  to  avoid  any  overlap  of  Images  between  neighboring  elements  on  the 
receiver  plane,  it  is  necessary  to  have  an  opaque  area  between  neighboring 
elements.  The  width  of  this  area  is  determined  by  the  resolution  of  each 
component  of  the  system.  Similarly,  when  each  light  source  is  not  an  ideal 
point  source,  the  size  of  the  light  sources  and  the  gap  between  them  are 
deterained  by  the  resolution  of  the  system  and  the  element  size  on  the 
encoding  plates.  Because  optical  wave  theory  can  only  provide  a  rough 
approximation  with  regard  to  the  data  (and  it  often  is  far  off  from  reality) , 
it  will  not  be  discussed  in  detail.  Experimentally,  the  spacing  between  two 
neighboring  elements  was  chosen  to  be  R  times  the  width  of  the  element  itself 
(when  approaching  the  resolution  limit  of  the  system,  R  >  1;  when  far  away 
from  the  resolution  limit,  R  «  1). 
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2.  Error  due  to  nonuniform  angular  light  distribution  and  ways  to  eliminate  It 

The  error  due  to  nonuniform  angular  distribution  of  a  light  source  Is  obvious. 
With  Gaussian  distribution,  It  can  be  expressed  as: 

J(0)  =  Joe-**” 

When  the  light  from  a  point  source  passes  through  the  matrix  multiplier  and 
reaches  A,  it  can  be  expressed  as  follows: 

•‘■AY  ® 

y  -  /  tg  «  /« 

After  penetrating  plate  A,  It  Is 


Upon  reaching  B,  It  becomes 

After  penetrating  B,  It  Is  ^  receiver  plane,  we  have 

Ixiyi  “ 

where  G  is  dependent  upon  the  diaphragm  of  the  encoding  plate. 

In  order  to  allow  this  distribution  to  meet  the  accuracy  requirement  we 
have 


f/D  -  [18a/ln(m)]i/2 

where  a  Is  a  parameter  which  defines  the  Gaussian  distribution,  5  Is  the 
accuracy  requirement ,  and  D  Is  the  maximum  width  of  the  encoding  plate ,  1 ,  e .  , 
D  -  max  (Dx,  Dy) .  In  practice,  different  light  sources  may  be  chosen  based  on 
the  data.  Of  course.  It  Is  only  feasible  based  on  availability. 

3.  Methods  to  eliminate  error 

Spatial  distribution  multiplication  error  can  be  eliminated  by  a  simple 
method.  It  involves  adding  encoding  plates  with  Ay  -  Cy  =  1  and  Sy  -  1  onto 
matrix  plates  A  and  C.  Use  the  receiver  to  measure  and  determine 
min(2>^^)/p^_^  -  Dy.  Place  an  encoding  plate  with  transmittance  Dy  immediately 
in  front  of  the  receiver.  This  is  a  way  to  eliminate  any  multiplication  error. 

Similarly,  additive  error  may  be  eliminated  electronically  at  the  receiver. 
Using  this  method,  error  can  be  minimized  to  a  very  low  level  (below  1  per¬ 
cent)  ,  because  most  errors  in  the  optical  system  can  be  eliminated  by  these 
two  techniques. 
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System  Adjustment,  Experimental  Procedure,  and  Results 

Using  the  system  shown  in  Figure  1,  experiments  were  done  to  multiply  two 
matrices  Aioxs^sxio  three  matrices  U^oxsi  ^5^5,  and  W5X10. 

1.  System  adjustment 

First,  align  the  center  of  each  lens  and  encoding  plates.  Adjust  each  lens  so 
that  it  is  perpendicular  to  the  center  line.  Then,  adjust  the  two  sets  of 
cylindrical  lenses  to  make  them  orthogonal  to  each  other.  Based  on  the 
dimension  and  accuracy  requirement  of  each  optical  element  and  encoding  plate, 
the  following  parameters  can  be  obtained: 

center  deviation  <  (d-D)/4; 

vertical  deviation  <  h/D/2 

where  d  is  the  actual  aperture  of  the  optical  element,  D  =  max(Dx,  Dy) ,  and  h 
is  the  minimum  width  of  an  encoding  plate  element. 

2.  Experimental  procedure 

In  the  experiment,  /  -  10  mm,  each  encoding  plate  element  is  approximately 
0.6  X  0.6  mm^,  and  the  gap  is  0.4  mm  or  0.3  mm  wide.  Every  cylindrical  lens  is 
mounted  on  a  two-dimens ionally  adjustable  rack.  Each  encoding  plate  is  clamped 
onto  a  three-dimens ionally  adjustable  rack.  Experimentally,  five  point  sources 
in  were  first  used  for  illumination.  After  using  encoding  plates  A,‘j  -  1  and 
C\j  -  1,  adjustments  were  made  to  obtain  uniform  results  on  the  receiver 
surface.  Then,  the  real  encoding  plates  were  inserted  to  conduct  the  tests. 
The  error  of  the  display  is  1.0  percent. 

3.  Experimental  results 

(1)  Multiplication  of  two  matrices; 


110  0  1' 
1110  0 
0  1110 
0  0  111 
10  0  11 
110  0  1 
1110  0 
0  1110 
0  0  111 
.10  0  11. 


1  0  0  0  0 
0  10  0  0 
0  0  10  0 
0  0  0  1  0 
0  0  0  0  1 


1100111001 

1110011100 

0111001110 

0011100111 

1001110011 


A  B 


c 
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■3  21123211  2' 
232112  3  2  11 
1232112321 
1123211232 
2112321123 
3211232112 
2321123211 
1232112321 
1123211232 
.2  11232112  3. 

(2)  Multiplication  of  three  matrices: 

1  1  0  0  O' 

1110  0 

0  1  1  1  0  r  1  r 

100011000110001 

0  0  111 

11000  1100011000 

11001  ^  ^  ^  ^  ^  011000  1  1  00 
00110  0011000110 

1110  0 

^  00011  0001100011 
0  1  1  1  0 

0  0  111 

10  0  11. 

A  B  C 

■3  101331013' 
4311343113 
3431134311 
1343113431 
1134311343 
ASC  3113431134 

4311343113 
3431134311 
1343113431 
1  13431134  3. 


Figure  2  [photo  not  reproduced]  shows  the  result  of  multiplying  two  matrices. 
The  picture  was  enlarged  (and  a  correction  plate  as  the  one  described  earlier 
was  inserted  in  the  experiment).  For  convenience,  each  element  is  a  plate  with 
either  1.0  or  0.0  transmittance.  Table  1  shows  the  optoelectronic  measurement 
results  for  the  situation  shown  in  Figure  2  [photo  not  reproduced] . 
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Table  1.  Results  of  Two-Matrix  Multiplication  (Detected  by  optoelectrical 
detector) 


8.01 

2.02 

0,98 

1.01 

1.98 

8.03 

2.04 

0.99 

1.03 

2.02 

1.97 

2.9S 

1.97 

1.02 

0.99 

1.99 

2.99 

2.01 

0.98 

0.99 

2.02 

8.03 

2.01 

1.01 

0.98 

2.08 

8.M 

2.03 

0.99 

1.03 

1.01 

2.04 

8.05 

1,97 

0.99 

1.01 

2.03 

3.01 

1.97 

1.98 

0.98 

1.01 

2.03 

2.97 

1.99 

0.99 

1.01 

1.99 

3.03 

8.06 

2.01 

0.99 

0.99 

2.01 

2.97 

1.96 

0.98 

1.01 

1.98 

2.02 

3.01 

2.03 

1.02 

1.01 

1.99 

2.96 

2.02 

0.98 

0.99 

0.99 

1.97 

2.95 

1.97 

0.98 

0.99 

2.03 

3.05 

2.03 

1.01 

1.00 

0.99 

2.01 

2.95 

2.04 

1.01 

1.02 

2.03 

3.03 

2.04 

1.99 

0.97 

1.01 

2.01 

3.02 

2.03 

1.01 

1.00 

2.01 

3.04 

Error-.. —2%  ETiorM»~1.2% 
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[Text]  Abstract 

Optical  multistability  and  optical  ternary  logic  gates  have  been  realized 
using  a  composite  ZnS  interference  bistable  device.  The  device  is  an  ideal 
element  for  studying  the  feasibility  of  optical  computing  because  of  its 
favorable  characteristics  such  as  ease  of  fabrication,  low  threshold  power, 
and  room-temperature  operation. 

Interference-filter-based  optical  bistable  devices  are  simple  to  fabricate  and 
easy  to  operate.  They  may  be  used  as  logic  elements  in  areas  such  as  optical 
communications,  optical  computing,  and  optical  information  processing. 
Interference  filter-based  optical  bistable  devices  have  been  successfully  used 
in  optical  full  adders,^  mode  recognition, ^  digital  conversion,®  and  high-speed 
switching.^  This  paper  reports  the  use  of  a  composite  interference  filter 
bistable  device  to  realize  optical  multistability  and  optical  symmetric 
ternary  logic  gates.  Compared  to  binary  logic  systems,  a  multi-valued  logic 
system  has  higher  operating  speed  and  information  storage  density.  Recently, 
Liu  Shutian  and  M.  Watanabe  investigated  optical  multi-valued  logic  devices 
such  as  the  Post  ternary  logic  gate  and  S3mimetric  ternary  logic  gate  using 
optoelectronic  mixers,®*®  respectively.  Their  work  has  considerable 
significance  in  improving  the  capability  of  optical  computing. 

I.  Optical  Multistabllity 

A  dielectric  interference  filter  contains  alternate  coatings  of  high- 
refractive  index  material  (such  as  ZnS  or  ZnSe)  and  low— refractive  index 
material  (such  as  NagAlFg  or  MgFj) .  Each  layer  is  one-quarter  wavelength  thick. 
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The  film  structure  is  (HL)'"(nH)  (LH)*".  (HL)"  and  (LH)"  are  considered  as  two 
highly  reflective  mirrors  on  both  sides  and  the  entire  filter  is  equivalent  to 
a  nonlinear  F-P  standard  device.  The  ZnS  interference  filter  is  intrinsically 
a  thermally  activated  optical  bistable  device.  The  bistability  is  caused  by  a 
shift  of  the  center  wavelength  of  the  filter  to  the  incident  wavelength  due  to 
heat.  The  threshold  power  to  produce  bistability  is  dependent  upon  the 
incident  angle.  The  larger  the  angle  of  incidence  is,  the  higher  the  threshold 
power  becomes . 


In  this  experiment,  two  interference 
filters  were  used  in  series  to  obtain 
optical  tristability.  The  experimental 
setup  is  shown  in  Figure  1 .  The  output  of 
the  first  interference  filter  is  used  as 
the  input  for  the  second  filter.  The  fil-  Figure  1.  Experimental  Setup 
ters  have  a  center  wavelength  of  513.1  nm  Ar"^:  argon  laser;  M:  modulator; 
and  513.8  nm,  a  FWHM  of  2.1  nm  and  1.7  nm,  BS:  beam  splitter; 
and  minimum  threshold  power  of  15  mW  and  lenses;  IFj,  IF2: 

11  mW,  respectively.  Figure  2  [photo  not  interference  filters; 


reproduced]  shows  a  tristability  curve  D^,  D2:  detectors 

obtained  by  varying  the  input  light  power. 


The  tristability  produced  is  in  reality 
the  superposition  of  two  bistable 
effects.  As  the  incident  light  intensity 
increases ,  the  output  light  intensity 
also  continues  to  rise.  When  the  input 
light  intensity,  P,-,  is  equal  to  that  of 


the  threshold,  IF^,  IF^  is  opened.  At 
this  time,  although  I F2  is  still  at  its 
lower  state,  the  output  light  intensity, 
Pq,  still  Jumps  considerably  because  of 
its  relatively  high  background  transmit¬ 
tance.  As  the  input  light  intensity  in¬ 
creases,  when  the  light  intensity  pass¬ 
ing  IFi  reaches  the  threshold  of  IF2,  IF2 
is  opened  and  the  output  li^t  intensity 
jumps  up  once  again.  This  is  how  tri¬ 
stability  is  generated.  The  shape  of  the 
curve  may  be  varied  by  changing  the 


Figure  3.  Truth  Tables  and  Input- 

Output  Relation  Curves  for 
"OR,"  "AND"  Gates 


incident  angles  of  IFi,  and  IF2. 


II.  Optical  Multi-Valued  Logic 


The  most  basic  operations  in  symmetric  ternary  logic  are  "AND"  and  "OR." 
Figure  3  shows  the  truth  tables  and  the  input-output  relation  curves  for  these 
two  logic  gates.  It  indicates  that  OR  and  AND  gates  can  be  realized  by  the 
same  optical  path.  It  only  requires  changing  the  thresholds  of  the  two 
bistable  devices.  Under  conditions  to  meet  the  logic  relation  for  the  AND 
gate,  it  is  possible  to  meet  the  logic  relation  for  the  OR  gate  by  raising  the 
threshold  of  the  first  bistable  device  and  lowering  that  of  the  second  device 
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at  the  same  time.  Experimental  results  of  the  AND  and  OR  gates  obtained  using 
the  optical  system  illustrated  in  Figure  1  are  shown  in  Figure  4  [photo  not 
reproduced] .  Pulses  of  different  amplitudes  correspond  to  different  input 
light  intensities.  Pq,  P^,  and  P2  are  the  discriminating  intensities  of  0,  1, 
and  2,  respectively. 

The  NOT  gate  [or  "inverter"]  is  also  a  common 
logic  gate.  It  may  also  be  realized  with  two 
filters ,  as  shown  in  Figure  5 .  The  light 
reflected  by  the  second  filter  is  the  output. 

Due  to  polarization  Pb,  when  the  input  is  0,  1, 
and  2,  the  output  is  2,  1,  and  0.  This  meets 
the  logic  relation  for  the  NOT  gate. 

III.  Discussion 

Interference  filters  are  used  to  study  the 
feasibility  of  optical  computing.  Later,  they 
can  be  replaced  by  high-speed  switches.  Because 
optical  systems  have  a  high  degree  of  parallel¬ 
ism,  high-processing-speed  optical  computing  ,  ^  „ 

systems  of  practical  value  can  be  achieved  even  5.  Principle  of  NOT 

with  medixam-speed  logic  elements  (1  /is-l  ns)  .  Operation 

The  use  of  medium-speed  logic  elements  is  also 

compatible  with  available  light  sources  and  signal  detectors.  As  progress  in 
interference  filter  bistable  devices  advances,  ZnSe  filters  operating  in  the 
infrared  have  been  successfully  developed.  Such  devices  can  be  used  with  laser 
diodes  to  substantially  reduce  the  volume  of  the  system. 
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[Text]  Abstract 

The  optical  parallel  cache  memory  array  and  random  access  array  are  critical 
parts  of  a  digital  optical  processor.  A  simple  optical  parallel  cache  memory 
array  Is  presented.  Furthermore,  random  access  to  the  optical  array  has  been 
achieved  experimentally. 

I .  Introduc t Ion 

Optical  mathematical  computation  Is  a  highly  parallel  spatial  optical  logic 
operation.^  It  allows  the  simultaneous  transfer  and  processing  of  a  vast  amount 
of  optical  Information  through  a  hlgh-denslty  parallel  optical  network  and 
two-dimensional  optical  switching  array. ^  From  the  algorithm  and  system 
architecture  of  digital  optical  computing,  the  two-dimensional  optical  logic 
array  Is  essential  to  such  a  parallel  digital  optical  processor.  The  optical 
parallel  cache  memory  array  and  parallel  random  access  memory  array  are  also 
vital  to  such  a  system.^  Because  digital  optical  computing  Is  highly  parallel. 
In  order  to  suit  high-speed  optical  processing,  the  data,  command,  and  program 
required  for  processing  are  first  stored  In  the  cache  memory  array  for  ready 
access.  In  optical  Information  processing,  there  Is  also  a  need  to  store 
Intermediate  results  In  random  access  memory  to  be  taken  out  In  the  next  step.^ 

There  are  numerous  studies  on  optical  bistable  logic  devices. In  this  work, 
a  random  access  memory  array  consisting  of  two-dimensional  optical  bistable 
devices  Is  presented.  Details  regarding  writing,  reading,  and  zero  setting  are 
described.  Parallel  optical  random  access  of  Information  was  experimentally 
achieved. 
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II.  Optical  Bistable  Memory 

In  digital  optical  computing,  an  optical  bistable  device  can  be  used  as  an 
optical  logic  switch,  as  well  as  optical  cache  memory  array. ^  As  an  optical 
logic  device,  its  bistability  curve  should  be  adjusted  to  approach  0,  as  shown 
In  Figure  1.  As  an  optical  memory  device.  It  requires  a  more  stable 
bistability  range  in  order  to  implement  basic  functions  such  as  random 
writing,  storage,  reading,  and  clearing. 


I  In  digital  optical  ■ 

Pout  computing,  the  binary 

J— ^  code  of  spatial  logic 

*  consists  of  a  series  ^  _ 

of  black  and  white  ~~  c 

b  patterns. A  bright  f 

- spot  has  a  certain  I  I 

Figure  1.  Output/Input  —7=  ^  ^ 

1  V-  power,  approach-  Figure  2.  Outpuc/Input 

logical  Switching  »  Characteristic 

of  discussion,  let  us  for  the  Memory 

assume  a  bright  spot  is  "1"  and  a  dark  spot  is 
"0."  When  a  bistable  device  is  used  as  a  memory 

device,  as  shown  in  Figure  2,  the  power  of  bias  light  Hi  is  first  adjusted  to 
point  b  in  the  bistable  zone.  When  the  input  signal  is  "0,"  point  b  belongs  to 
the  lower  state  and  the  output  signal  is  "0."  When  the  input  signal  is  "1,"  it 
jumps  from  point  b  to  point  c  which  is  at  a  higher  state  and  the  output  is 
"1."  After  the  signal  is  removed,  the  output  still  remains  at  point  b' .  There¬ 
fore,  the  original  optical  signal  "1"  is  stored  in  memory.  If  this  information 
needs  to  be  randomly  accessible,  an  additional  optical  bistable  switch  is 
required.  Its  characteristic  curve  is  shown  in  Figure  1.  The  bias  light  for 
this  device  is  set  at  point  b  before  the  bistable  zone.  If  one  needs  to  take 
out  the  information  at  a  certain  point  of  the  memory,  the  bias  light  must  be 
applied  to  the  corresponding  point  on  the  switch  to  read  it  out. 


Figure  1.  Output/Input 

Characteristic  for 
Logical  Switching 

assume  a  bright  spot  is 


Figure  3 .  The  Experimental  Setup 
of  4  X  2  Array  for  Parallel 
Optical  Cache  Memory  Array 


If  the  next  piece  of  information  needs  to 
be  stored  at  the  same  location,  then  the 
data  stored  there  must  be  cleared.  This  can 
be  done  by  cutting  off  the  bias  light  on 
the  memory  and  the  corresponding  bias  light 
H2  on  the  switch  to  set  0. 

III.  Optical  Parallel  Cache  Memory  Array 
Experiment 

Figure  3  shows  the  experimental  setup  for 
the  optical  system.  The  optical  bistable 
device  is  the  normal  transmission-type 
optical  bistable  ZnS  device®*®  developed  in 
our  laboratory.  Its  bistability  character¬ 
istic  curve  is  shown  in  Figure  4. 
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The  X  and  y  axes  are  the  input  and  output  power,  respec- 
tively.  In  the  optical  path,  the  planar  thin-film  beam 
splitter  splits  the  beam  into  a  4  x  4  array.  The  SSR 
composite  prism  and  the  4x4  microlens  array  were  all 
developed  in  our  laboratory. 

The  SSR  composite  lens  can  divide  the  4x4  beam  array 
into  two  arrays  of  identical  intensity.  Moreover,  the 
two  beam  arrays  can  be  displaced  relative  to  each  other 
and  make  beams  from  different  points  in  space  coincide 
with  one  another  in  a  rigorous  manner.  In  this  experi¬ 
ment,  the  first  and  third  column  and  the  second  and 
fourth  column  were  made  to  overlap  by  the  SSR  prism  to 
form  a  4  X  2  array,  as  shown  in  Figure  3.  Using  the  half 
wave  plate  HWPi,  the  power  of  each  beam  in  column  P  is 
adjusted  to  9  mW.  P^  and  P2  are  used  as  the  bias  light 

Hi  for  the  memory.  In  the  S  light  array,  the  S3S4  array  is  the  signal  light 
source.  Each  beam  has  a  power  of  3  mW.  Arrays  S3S4  and  P1P2  are  overlapped  and 
passed  through  the  microlens  to  enter  the  normal  transmission  optical  bistable 
device  OBDi  at  right  incident  angle.  Array  P3P4  is  used  as  the  bias  light  H2  for 
the  second  optical  bistable  device  OBD^.  Its  optical  power  may  be  adjusted  by 
HWPi. 


-m 


Figure  4.  Bistable 
Characteristic  of 
Normal  Transmission 
ZnS  IF  Device 


:4: 


(a1  (b)  fc)  (d) 

Figure  5.  The  Process  of  Writing  to,  Reading  From, 
and  Clearing  the  Memory  in  Optical  Cache 
Memory  Array 


In  the  experiment,  a 
moving  screen  M^  was  used 
to  block  off  light  beams 
in  columns  3  and  4. 

Therefore ,  P1P2  is  the 
only  bias  light  shining 
on  OBDi.  If  the  bias 
light  is  modulated  sinu¬ 
soidally  from  0  to  9  mW, 
the  bistable  curve  of 
the  4x2  array  output 
of  OBDi  is  as  shown  in 

Figure  5(a).  From  the  photograph,  when  the  bias  light  power  is  at  9  mW  for 
every  beam,  every  point  of  the  4x2  array  is  at  its  low  state.  When  the  input 
light  signal  is  "0,"  the  OBDi  output  signal  is  "0."  When  M2  is  used  to  block 
off  the  light  array  P3P4  and  Mi  is  removed  to  let  the  light  array  838^  with 
optical  signal  "1"  through,  every  point  of  the  4x2  light  input  array  for  OBDi 
simultaneously  jumps  from  "0"  to  "1,"  as  shown  in  Figure  5(b).  If  the  incident 
light  (bias  and  signal)  is  sinusoidally  modulated,  i.e.,  the  incident  light 
power  is  varied  from  0  to  12  mW  sinusoidally,  the  bistable  curve  for  every 
point  of  the  4x2  input  array  for  OBDi  shown  in  Figure  5(c).  The  photo¬ 
graph  indicates  that  when  the  bias  is  9  mW  and  the  signal  is  3  mW,  all  points 
in  the  4x2  OBDi  output  array  are  at  their  high  states  and  the  output  is  "1." 
Using  Ml  to  block  off  8384,  the  4x2  output  array  of  OBDi  is  still  at  the  high 
state,  i.e.,  a  bright  spot  "1"  as  shown  in  Figure  5(d).  This  illustrates  that 
a  4  X  2  input  array  of  "1"  is  written  into  and  stored  in  memory.  However, 
there  is  no  light  signal  output  from  OBD2.  When  the  optically  stored  informa¬ 
tion  is  needed,  screen  M2  shown  in  Figure  3  is  removed  to  allow  the  bias  H2  to 
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shine  on  the  4x2  array  in  OBD2  in  the  form  of  array  P3P4.  The  stored  data  can 
then  be  read  immediately  because  every  point  in  OBD2  jumps  from  a  low  state  to 
a  high  state  and  the  output  optical  signal  is  "1.” 


If  we  want  to  store  the  information  from  the 
next  interval  in  the  same  point  of  memory, 
the  data  stored  must  be  cleared  first.  Before 
storing  new  data,  the  bias  light  must  be 
removed  from  the  point  to  be  cleared.  After¬ 
ward,  bias  light  can  be  restored  to  begin 
storing  new  data.  The  writing,  reading,  and 
clearing  of  data  are  shown  in  Figure  6. 


^set  zero 

^writing  pulse 
output  of 


n 


reading  pulse 


output  of  OBDo 


From  this  experiment,  we  know  that  the  system 
can  write,  store,  read,  and  clear  data.  In 
optical  computing,  it  can  serve  as  a  parallel 
cache  memory  array,  as  well  as  a  parallel 
random  access  memory. 


Figure  6.  The  Process  of 

Writing  Into  and  Reading 
From  the  Optical  Cache 
Memory 


Conclusions 


The  experimental  optical  parallel  cache  memory  array  consists  of  two  uniform 
normal  bistable  ZnS  devices.  The  optical  system  includes  a  high-efficiency 
planar  thin-film  beam  splitter  to  provide  a  uniform  beam  array.  A  multipurpose 
SSR  prism  ensures  the  stability  of  the  optical  system.  This  simple  optical 
system  can  perform  basic  parallel  functions  such  as  write,  store,  read,  and 
clear  in  a  random  fashion.  The  simultaneous  access  of  an  optical  signal  array 
is  primarily  determined  by  the  homogeneity  of  the  beam  array  and  the  uniform¬ 
ity  of  the  bistable  device.  Various  parameters  of  the  optical  cache  memory 
array  depend  upon  the  upward  or  downward  transition  time.  The  reliability  of 
optical  storage  depends  upon  the  stability  of  the  bistable  zone.  In  digital 
optical  computing,  this  system  can  be  used  as  a  parallel  cache  memory  array  of 
a  parallel  random  access  memory  array. 
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[Text]  Abstract 

Matrix-vector  and  matrix-matrix  multiplication  have  been  achieved  using  four- 
wave  mixing  in  LiNbOstFe  crystal.  Light  intensities  representing  matrix 
elements  are  below  100  /iW/mm^  and  the  response  time  is  of  the  order  of  seconds. 
The  photorefractlve  effects  affecting  the  accuracy  of  the  operation  are  also 
discussed. 

I .  Introduction 

Because  all-optical  computation  based  on  photorefractlve  nonlinear  optics  is 
parallel  and  can  handle  data  over  a  wide  frequency  band,  it  is  a  very  attrac¬ 
tive  subject.  Yeh  and  Chiou^  first  accomplished  pixel-pixel  multiplication 
using  four-wave  frequency  mixing  in  BaTiOs  in  1987.  Then,  these  pixel  products 
were  summed  with  a  cylindrical  lens  to  complete  the  matrix-vector  multiplica¬ 
tion  experiment.  This  kind  of  vector-matrix  multiplication  can  be  extended  to 
matrix-matrix  multiplication  by  using  composite  light. ^  Matrix-matrix  multipli¬ 
cation  can  also  be  done  by  finding  the  spatial  convolution  of  four-wave 
frequency  mixing.  Chiou  and  Nguyen®  implemented  this  matrix-matrix  multiplica¬ 
tion  scheme  with  BaTlOs.  As  far  as  optical  computing  is  concerned,  the  light 
intensity  response  of  the  photorefractlve  crystal  must  be  sufficiently  low  so 
that  a  large  number  of  light  beams  can  be  produced  by  a  laser  to  complete  the 
multiplication  of  higher-order  matrices.  Therefore,  finding  a  crystal  that  has 
a  low  response  to  light  intensity  is  critical.  It  was  experimentally  found 
that  the  optical  response  of  LiNb03  with  0.1  wt%  (weight  percent)  Fe  is  below 
100  /iW/mm®.  Furthermore,  the  photorefractlve  effects,  i.e.,  light-induced 
scattering*  and  self-defocusing,®  have  substantial  impact  on  accuracy. 
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II.  Matrix-Vector  Multiplication 

1.  Multiplication  by  Four-Wave  Mixing  in  Photorefractlve  Crystal 

The  four-wave  mixing  process  in  a  photorefractlve  crystal  can  be  described  by 
a  steady-state  wave-coupling  equation.  As  shown  in  Figure  1,  four  beams  of 
light  with  frequency  w  and  Identical  polarization  intersect  in  the 
photorefractlve  crystal.  Let  the  photoelectric  field  of  these  four  beams  be: 


Ej(r,t)  “  Aj(r)exp[i(kj -r-wt) ]  +  c.c., 


(1) 


where  j  =1,  2,  3,  4.  Furthermore,  the  wave  vectors 
satisfy  these  relationships:  ki  -  — k2,  k3  —  — k^. 
These  four  beams  induce  four  refractivity  gratings 
in  the  crystal.  Two  of  them,  i.e.,  transmission 
grating  and  reflection  grating,  are  the  medium  of 
four-wave  mixing.  The  wave-coupling  equations 
resulting  are  as  follows®: 


(1)  The  four-wave  mixing  coupling  equations 
primarily  due  to  the  transmission  grating  are: 


=  ^(AiA;+A2*A3)A2+OA3*, 

dz  «X0 


=  ^(A3a;+a2*a3)a;-«a;, 
az  Jo 

iMiij  exp(-i4>j)  ^  ^  I, 

2c  cos  0  ~  ^ 


(2) 


^2  z=l  2=0 


Figure  1.  Diagram  of  the 
Four— Wave  Mixing 


where  nj  is  the  refractivity  grating  amplitude,  is  spatial  phase 
difference  between  the  grating  and  the  interference  fringe,  c  is  the  speed  of 
light,  and  6  is  the  angle  between  the  light  velocities. 


(2)  The  four— wave-mixing  coupling  equations  primarily  due  to  the  reflection 
grating  are: 


dAi 

dz 

dA^ 

dz 

dAl 

dz 

dA^ 


dz 

Y 


=  -3L(A3A3*+A2*A.)A3-«A3, 
■^0 

=  ■^{\At*A2A^)A^*nAi, 

1q 

=  -^UtA;+A;A^)Al*aA;. 

Iq 

=  ^(AiA3*+A2*A.)A2-«A., 
•^0 

i<anjj  exp(i<|>jj) 

2c  cos  6  ’  j 


(3) 


22 


In  LiNbOsrFe  crystal,  the  primary  mechanism  for  electron  transport  is 
diffusion.  Therefore,  (f>i  -  -  (,n/2)  Because  7  is  a  real  number  and  when 
there  is  almost  no  absorption  (a  -  0) ,  these  two  sets  of  equations  can  be 
rigorously  solved.® 


When  equation  (2),  i.e.,  transmission  grating,  plays  a  dominant  role,  the 
conjugate  reflectivity  is 


■^3(0)  ^  _  I _ -2cT 

A^iO)  '  I  Ar+(A»+4|c|2)i/2 


(4) 


A^{t)  =  0,  J1+I4  ®  const  =  Ji(0)+J4(0),  J2+X3  =  const  =  JjCt), 
A  =  {J2+X3)  -  -  A(o)  -  •r4(o)» 

J3  +  J4  =  Ji(0)  +  J2(l)  +  lt(0)  . 

Q=  (A*+4  p  =  (yQ/2If,) ,  r=  tanhpl, 

|c|^  is  determined  by  the  following  equation: 


[1C|2  -  Ii(0)l2(i)]  |AT  +  (a2+4|C|2)1/2|2  +  4|Cl2|T|2l^(0)l2(je) 


+  2|C|2i^(0)  (a2+4|C|2)1/2(t+t*)  -  0. 

Figure  2.  Dependence  of  the  Phase- 
Conjugate  Reflectivity  on 
the  Pumping  Intensities.  The 
theoretical  fitting  curves 
are  obtained  for  7^  --1, 

I4  -  1.0  mW.  The  numerical 
values  above  each  curve  cor¬ 
respond  to  the  intensities 
of  Ii  (xlOO  pW). 


(6) 


For  I4  =  1.0  mW,  Ij,  I2  <  100  pW,  theoretical  fitting  curves  for  R  versus  li(0) 
and  l2(-?)  have  been  computed  and  are  shown  in  Figure  2.  When  I4  »  I^,  I2  and 
I3,  we  have 

A  i  -l4(0),  Q  i  14(0),  M  i  (7/2),Io  -  Ia(0),  (7) 

Equation  (6)  can  then  be  simplified  as 

|C|2  -  [(1-T)/(1+T)]2li(0)l2(i).  (8) 

Furthermore ,  from  equation  (4) ,  we  get 

^  ~  f  1  J3 (0)  J2  »)  «Ji  (0)  J2  (t)  .  (9) 

^  J4(0)exp(Y«)  J  1212 

When  equation  (3),  i.e.,  reflection  grating,  is  dominant,  its  conjugate 

reflectivity  is 
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tanhM  1^0 !>»«)]  “  v^v  . 


_ Ji(0)  J4(0)  [1-exp (-Yt)]^ _ 

{Ii  (0)  +  [Jj  («)  +  J4  (0)  ]  exp(-2l)  }*+Ji (0)  I4 (0)  [l-exp(-Y«)  ]  *  ' 


For  LiNbOsrFe,  |7i|  ~  1.  Thus, 


tanh2[|go||i(l)] 


^  Ji(0)  [1-exp (-yt)]' 
14(0)  exp  (-2y<) 


1-exp  (-Y<) 

J4(0)exp(-Y«) 


From  equations  (9)  and  (11),  we  can  see  that  either  transmlsslon-gratlng-  or 
reflectlon-gratlng-lnduced  four-wave  mixing  can  accomplish  multiplication  as 
long  as  the  probe  ratio  q  -  {14(0)  [Ii(0)+l2(i)  ] }  »  1.  When  yi  <  0,  the 
transmission  grating  plays  the  major  role.  However,  when  yi  >  0,  the 
reflection  grating  dominates.  The  use  of  yi  Is  determined  by  the  crystal  axis 
and  polarization  of  the  coupled  light  (l.e.,  the  effective  photoelectric 
coefficient  yeff®’^)  • 

2.  Matrix-Vector  Multiplication  by  Image -Beam/Image -Beam  Multiplication 


photoreftaetive 

crystal 


probe  beaa 
(plane  wave) 
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Figure  3.  Conceptual  Diagram  of 
Matrix-Vector  Multiplication 


Figure  4.  Experimental  Setup  for 

Matrix— Vector  Multiplication 
Ml,  M2,  M3  are  mirrors;  Li,  L2,  L3  are 
cylindrical  lenses;  BS,  BSi,  BS2  are 
beam  splitters.  CL  Is  the  collimator. 


Figure  3  shows  the  principle  of  matrix-vector  multiplication.  The  matrix  and 
vector  elements  are  coded  by  two  aperture  masks  with  light  passing  in  opposite 
directions .  A  collimated  Intense  beam  of  light  Is  used  as  the  probe .  Experi¬ 
mentally,  a  1.5-mm-thick  LiNb03:Fe  crystal  is  used  for  the  image-beam/image- 
beam  multiplication.  These  products  are  "read"  by  the  probe  light,  which  is 
focused  by  a  cylindrical  lens.  The  summation  is  done  on  the  focal  plane  to 
complete  the  matrix-vector  multiplication  process.  The  experimental  setup  is 
shown  in  Figure  4.  Figure  5  shows  the  result  of  the  product  of  a  2  x  2  matrix 
and  a  2  X  1  vector.  For  simplicity,  a  thick  beam  was  used  in  the  experiment, 
instead  of  using  cylindrical  lenses  Li  and  L2  to  expand  the  beam. 


lij]  l|]  over 
ITi  U2  ^1^2 


Figure  5.  The  Operation  of 
Matrix-Vector  Multiplication 

III.  Matrix-Matrix  Multiplication 


Figure  6.  Conceptual  Diagram 
of  Matrix  Multiplication 
Using  Convolution 


1.  Convolution  and  Correlation  Using  Four -Wave  Mixing 

White  and  Yariv  proved  that  the  correlation  and  spatial  convolution  of  two 
patterns  could  be  done  by  four-wave  mixing.^®  The  crystal  used  in  that 
experiment  is  Bli2Si02o  and  its  principle  is  shown  in  Figure  6.  After  Fourier 
transformation,  the  object  plane  input  amplitude  Ui(x,y),  U2(x,y),  and  UsCx.y) 
become  u^,  U2,  and  U3  (U3  »  u^,  U2)  on  the  confocal  plane.  These  three  beams 
are  coupled  in  the  crystal  to  produce  four-wave  mixing.  The  conjugate  of  U3, 
U4,  is  converted  to  U^Cx.y)  by  inverse  Fourier  transformation  using  lens  Ig. 
It  can  be  expressed  as  follows^”: 

U4(x,y)  -  ^U3(x,y)  ★  [Ui(x,y)  •U2(x,y)],  (12) 

where  is  a  transformation  constant,  ★  and  •  are  correlation  and  convolution, 
respectively. 

When  Ui  and  U2  are  the  matrix  masks  and  U3  is  an  aperture  (5  function),  we  have 


C(x,y)  «  U^{x,y)  •  U2^x,y)  =  U2(x-5,y-ti)d{cfn 

-w 

a.i'i' 


U4  (x,  y)  =  (x,  y)  ★  [ Ui  (x,  y)  (x, y)  ] 

m 

=  ffJ‘[/2(i,r])C(x+^,y+r))c^dH 

-m 

m 

=  C(x+li,y+i)}dlidn 

=  i|fC(x,y) , 


(13) 


(14) 
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where  (5]^  t^ij)  is  the  coordinate  of  element 
ij  on  matrix  mask  Up.  Equation  (14)  shows 
that  after  passing  through  lens  L3,  the 
output  amplitude  of  the  conjugate  wave 
produced  by  four-wave  mixing  on  the  focal 
plane  is  the  convolution  of  Ui  and  U2. 

If  Ui  and  U2  are  the  matrices  shown  in  Figure 
7(a)  and  (b) ,  then  their  convolution  is  as 
shown  in  Figure  7(c).  The  center  column  of 
the  convolution  corresponds  to  the  elements 
of  U1U2. 


2.  Matrix-Matrix  Multiplication  by  Spatial 
Convolution 
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Figure  7.  The  Convolution  of 
15 1  and  U2 


Figure  8  shows  the  experimental  setup.  The  horizontal  distance  of  holes  on  Ui 
and  U2  is  2a  =  1.5  mm,  and  the  vertical  distances  are  2b  =  3  mm  and  2c  =  1.5 
mm.  Figure  9(a)-(c)  shows  the  experimental  results  of  matrix  multiplication. 
[The  cross  terms  of  the  convolution  (left  and  right  column  in  Figure  7(c))  are 
not  shown  because  they  have  no  contribution  to  the  matrix  multiplication.] 


Figure  8.  Experimental  Setup  of  Matrix- 
Matrix  Multiplication 
Ml,  M2,  M3  are  mirrors;  Li,  L2,  L3  are 
spherical  lenses  with  equal  focuses; 
Ui,  U2,  U3  are  matrix  masks 


Figure  9.  The  Operation  of 
Matrix-Matrix  Multiplication 


IV,  Experimental  Results  and  Discussion 


Major  factors  affecting  accuracy  are  self-defocus ing  and  light-induced 
scattering.  Self-defocusing  diffuses  the  output  light,  as  shown  in  Figure  5. 
For  long  exposure,  the  light  spots  become  undistinguishable .  Self-defocusing 
is  caused  by  an  index  change  induced  by  a  Gaussian  beam,  as  shown  in  Figure 
10(b).  If  a  light  beam  is  uniform  and  weak,  then  the  index  of  refraction 
change  is  flat  and  small.  Consequently,  the  self-defocusing  effect  is  negligi¬ 
ble.  Light-induced  scattering  causes  the  conjugate  wave  to  lose  some  energy. 
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As  shown  in  Figure  10(a),  multiplication  is 
accurate  in  the  stable  region  marked  in  the 
figure.  If  scattering  is  strong  enough,  this 
stable  region  will  disappear.  Light- induced 
scattering  is  dependent  upon  the  diameter 
and  intensity  of  the  beam.  Zhang  Guangyin^^ 
pointed  out  that  there  is  a  threshold  dia¬ 
meter  for  this  type  of  scattering.  There¬ 
fore,  a  low-intensity  and  small-diameter 
beam  can  significantly  reduce  the  effect  of 
light- induced  scattering  and  self-defocus ing. 
In  this  experiment,  the  image  beam  intensity 
is  less  than  100  /iW/mm^  and  the  diameter  is 
less  than  0.1  mm.  A  stable  area  of  more  than 
two  seconds  was  obtained  experimentally.  In 
addition,  in  matrix-matrix  multiplication  by 
convolution,  long-time  exposure  will  lead  to 
noise  diffraction  spots,  as  shown  in  Figure  9. 
These  noise  diffraction  spots  are  higher-order 
diffraction  points  of  the  transmission  or  re¬ 
flection  grating.  Hence,  whether  for  matrix- 
vector  or  matrix-matrix  multiplication,  there 
is  a  definite  time  within  which  the  result  is 
accurate . 


Figure  10.  Influences  of  Self- 
Defocusing  and  Light-Induced 
Scattering 

(a)  The  energy  loss  caused  by 
scattering. 

(b)  The  refraction  index 
change  caused  by  laser 
beam.  1.  Index  change 
caused  by  Gaussian  beam. 
2.  Index  change  caused  by 
weak  uniform  beam. 


V.  Conclusions 

All— optical  computation  using  photorefractive  nonlinear  optics  has  a  wide 
range  of  future  applications.  In  order  to  make  it  practical,  the  light 
intensity  response  of  photorefractive  materials  must  be  significantly  reduced. 
Experimentally,  it  was  determined  that  LiNb03  crystal  heavily  doped  with  Fe 
(0.1  wt%)  has  a  great  deal  of  potential  in  this  field. 

The  authors  wish  to  thank  Associate  Professor  Zhang  Chunping  [1728  2504  1627] 
and  Instructors  Li  Laobing  [2621  5071  0365]  and  Cui  Yufen  [1805  3768  5358]  for 
their  assistance  in  the  experimental  work. 
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[Text]  Abstract 

An  IPA  (Inter-Pattern  Association)  memory  based  on  Fresnel  hologram  is 
presented.  It  has  higher  recognition  capability  for  nonindependent  storage 
modes.  Using  Fresnel  holograms,  an  optical  IPA  memory  with  large  number  of 
neural  elements  can  be  implemented. 

I.  Introduction 

Neural  networks  have  been  an  active  area  of  research  in  recent  years .  Associa¬ 
tive  memory  is  the  most  critical  component  in  a  neural  network.  In  most 
optical  neural  networks,  Hopfield's  model  is  used  to  form  associative  memory 
However,  Hopfield's  model  does  not  consider  the  relation  between  various  modes 
in  storage.  It  only  stresses  the  inter-pattern  association  of  elements. 
Hopfield's  model  is  valid  when  storage  modes  are  independent  and  the  ON  and 
OFF  probability  for  each  element  is  equal.  Such  a  system  is  a  poor  model  for 
dependent  modes.  When  it  is  necessary  to  distinguish  similar  modes,  the 
difference  between  them  becomes  very  important.  Therefore,  when  we  develop  an 
IPA  memory,  the  similarities  and  dissimilarities  of  storage  modes  must  be 
considered  to  make  them  function  differently  in  the  association  process. 

In  1990,  Taiwei  Lu^  proposed  using  an  IPA  model  to  fabricate  an  associative 
memory.  An  optical  neural  network,  consisting  of  a  monitor,  lens  array,  spa¬ 
tial  light  modulator,  and  computer,  has  been  constructed.  The  IPA  associative 
memory  matrix  and  the  effect  of  the  lens  array  on  each  pixel  are  shown  on  the 
monitor.  The  results  indicate  that  the  address-finding  capability  of  the  IPA 
neural  network  is  substantially  higher  than  that  of  the  Hopfield's  model. 
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Especially  when  the  Hanming  distance  of  the  storage  modes  is  small,  it  is 
capable  of  finding  the  address  correctly.  Thus,  the  IPA  neural  network  has 
more  stable  capability  than  that  of  the  Hopfield's  model.  However,  because  of 
the  resolution  of  the  monitor,  the  number  of  pixels  stored  in  the  system 
cannot  be  very  large. 

A  hologram  associative  memory  is  particularly  suitable  for  the  storage  of  two- 
dimensional  pictures  with  multiple  pixels,  A  theoretical  analysis  of  the 
characteristics  of  Fresnel  holograms  is  presented  in  this  paper.  An  IPA 
associative  memory  was  fabricated  on  the  basis  of  this  analysis.  Furthermore, 
a  large-pixel  IPA  optical  associative  memory  has  to  be  implemented. 

II.  Theoretical  Analysis 

The  IPA  model  decomposes  every  mode  into  several  parts  based  on  how  one  mode 
coincides  with  other  modes  in  the  set  of  storage  modes.  Each  part  has  a 
different  function  in  the  associative  process.  A  unique  part  of  the  mode  has 
a  more  important  function  than  other  common  parts  in  association. 

For  example,  in  storing  three  modes,  when  some  part  of  the  input  address¬ 
finding  mode  is  ON,  if  the  corresponding  parts  of  the  three  modes  are  all  ON, 
then  the  input  mode  may  be  any  one  of  the  three  modes.  There  is  no  way  to  tell 
which  one;  we  can  only  say  there  is  an  input  mode.  Hence,  this  part  only 
stimulates  the  common  parts .  If  the  corresponding  parts  of  two  modes  are  ON 
and  that  of  the  third  mode  is  OFF,  then  the  input  mode  is  one  of  the  former 
two  and  it  is  not  the  latter.  However,  it  is  not  possible  to  distinguish 
between  the  two.  Then,  this  part  stimulates  the  common  part  of  the  former  two, 
but  suppresses  a  unique  part  of  the  latter.  If  only  the  corresponding  part  of 
one  mode  is  ON  and  those  of  the  other  two  are  OFF,  then  the  input  must  be  the 
former.  This  part  stimulates  all  the  parts  of  this  mode  and  suppresses  some 
part  of  the  other  two  modes.  Taiwei  Lu  proposed  a  triple  gray-scale  associa¬ 
tive  matrix.  The  elements  are  —1,  0,  and  1.  In  an  optical  neural  network,  it 
is  easy  to  implement  stimulation  and  suppression.  This  kind  of  associative 
matrix  is  relatively  easy  to  materialize  optically.  Nevertheless,  the  function 
of  each  part  in  the  associative  process  in  an  IPA  memory  is  proportional  to 
the  number  of  pixels  contained  in  this  part.  When  the  Hamming  distance  of  two 
storage  modes  is  small,  there  is  a  substantial  overlap  between  modes.  Their 
unique  parts  play  a  limited  role  in  the  associative  process.  It  is  more 
difficult  to  resolve  these  modes.  We  believe  that  we  should  add  a  gray  scale 
in  the  associative  matrix  to  balance  different  parts.  Since  the  system  chooses 
a  threshold  for  output,  the  gray-scale  requirement  is  not  very  stringent, 
which  is  easier  to  realize  optically. 

The  following  is  a  theoretical  analysis  of  IPA  memory  with  Fresnel  hologram. 
The  optical  path  used  to  record  the  Fresnel  hologram  is  shown  in  Figure  1 .  P2 
is  the  imaging  plane  of  plane  Pj  through  lens  Li.  P3  is  the  imaging  plane  of 
plane  P4  through  lens  L2.  During  recording,  a  transparent  plate  with  a  trans¬ 
missivity  of  a(p,q)  is  placed  on  plane  P2  to  come  in  contact  with  the  holo¬ 
graphic  mask.  It  is  illiaminated  by  another  focused  light.  Another  transparent 
plate  with  a  transmissivity  of  b(x,y)  is  placed  on  plane  P4.  A  diffuse  scatter¬ 
ing  screen  is  placed  in  front  of  it. 
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Let  us  assume  that  the  diffraction  field  of  b(x,y)  on  plane  P2  is  B(p,q).  In 
addition,  the  distribution  of  this  field  on  plane  P3  is  noted  to  be  D(x' ,  y') 
bi/2(x»  ^  y').  Then,  we  have 

Bip.q)  =  Jjoix', y')exp[-ik{d*p  ain  Q)] 

•exp  [iJc{  (p  cos  0-x')  *+  (g-yO  /2d]  dx'fy' . 


Figure  1.  Fabrication  of  Memory 
Matrix  for  an  IPA  Holo¬ 
graphic  Associative  Memory 


The  exposure  on  the  recording  medium  is 
proportional  to 


Figure  2.  IPA  Holographic 
Associative  Memory 


1  {exp[ik(p2+q2)/2/]  +  B(p-q)  }a(p,q)  1^.  (2) 

Here,  d  is  the  distance  between  the  holographic  mask  and  plane  P3,  K  is  the 
wave  number  of  the  laser,  /  is  the  focal  length  of  the  lens,  9  is  the  angle 
between  the  optical  axes  of  the  two  lenses,  and  D(x' ,  y')  is  the  diffuse 
scattering  screen.  Masks  representing  modes  a,-  and  their  corresponding 
stimulating  and  suppressing  parts  b|  are  placed  in  positions  a  and  b, 
respectively,  prior  to  exposure.  After  processing,  the  amplitude  of  the 
transmissivity  of  the  hologram  is  proportional  to  the  amount  of  exposure. 
Terms  contributing  to  first-order  holographic  diffraction  include: 


t{p,q)  =  ^  {e:qp[-ik(p2+g2)/2/]B^(p,g)}a*(p,g)  .  (3) 

Figure  2  shows  the  optical  path  of  a  holographic  IPA  memory.  After  placing  the 
recorded  Fresnel  hologram  on  plane  P2,  when  the  input  object  mask  a'(x,y)  is 
placed  on  plane  Pj  and  is  illuminated  by  a  collimated  beam,  the  transmission 
holographic  optical  field  is 

exp[-ik(p2+q2)/2/]a' (p,q) -tCp.q),  (4) 


On  output  plane  P3,  the  optical  field  can  be  expressed  as: 


?//*  i  (p,  g)  a'(p,  g)  (p, g) exp  [-iic(d+p  sin  0)  ] 
•  exp  [iJ:{  (p  cos  0-xO  *  +  (g-yO  /2d]  •  dpdq. 


(5) 
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Let  F|(x,y)  be  the  Fourier  transform  of  a,(p,q) -a' (p,q) .  Then,  the  output 
optical  field  can  be  written  as: 

{x-x' ,y-y')Di(x' ,y')bY^(x',y') 
exp  [iJc(jc*-jc'*+y®-y'*)/2ef]  •  c&'dy'. 

Because  of  the  presence  of  the  diffuse  scattering  screen  DjCx' ,  y'),  the 
intensity  distribution  can  be  expressed  as: 

't)ix-x',  y-y')bi^(x' ty'^dx'dy' 

=  ^  {a_,a'}i?i(x,y)  . 

The  bracket  {afa')  represents  the  inner  product  of  a  and  a'.  Here,  the  diffused 
light  spot  on  the  output  plane  is  neglected.  In  the  associative  process,  inner 
products  of  the  input  object  and  different  parts  are  obtained  to  determine  the 
corresponding  stimulation  and  suppression  distribution.  The  weighted  associa¬ 
tion  is  the  magnitude  of  the  inner  product.  The  magnitude  of  the  associative 
matrix  elements  can  be  implemented  by  adjusting  the  diffraction  efficiency  of 
the  holographic  mask.  In  this  work,  different  diffraction  efficiencies  were 
obtained  by  adjusting  the  exposure  time  to  record  the  holographic  mask. 

III.  Optical  Implementation 

Based  on  the  above  theoretical  analysis,  an  IPA  Fresnel  holographic  memory  has 
been  realized  optically.  Experimentally,  three  24  x  24  pixel  English  letters, 
X,  Y,  and  Z,  were  chosen  as  the  patterns  to  be  stored  in  binary  form.  Based  on 
this  concept,  pixels  can_be  divided  into  seven  categories:  (XnY)uZ,  (YnZ)uX, 
(ZnX)UY,  (XnY)nZ,  (Ynz)nx,  (Znx)n7,  (XnYnZ)nf,  where  n,  u,  and  —  represent 
logic  operations  AND,  OR,  and  NOT,  respectively.  ^  represents  the  set  of 
pixels  none  of  the  three  patterns  uses.  The  matrix  elements  with  suppressing 
effect  are  negative  in  value.  This  work  also  adopts  a  commonly  used  scheme  to 
provide  a  bias  for  every  matrix  element  so  that  the  minimum  occurs  at  zero 
transmissivity  in  order  to  compute  negative  values  optically. 

In  this  experiment,  a  30  mW  He-Ne  laser  was  used  as  the  light  source.  An  HP- 
366-1  holographic  mask  was  used  to  record  multiple  exposures.  Since  there  is 
no  overlap  between  the  light— transmitting  pixels  in  every  category,  the 
maximum  number  of  exposures  at  any  point  of  the  holographic  mask  is  limited  to 
one  to  ensure  the  diffraction  efficiency  of  the  mask.  A  two-dimensional  CCD 
detector  array  is  used  to  receive  the  output  picture.  The  signal  is  input  into 
a  computer  and  the  computer  selects  a  threshold  for  the  output  picture.  The 
associative  result  is  then  shown  on  the  monitor.  Figure  3(a)  shows  the  picture 
stored  in  the  system.  Figure  3(b)  is  a  part  of  the  input  stored  by  content- 
address  finding.  Figure  3(c)  is  the  associative  result  with  this  kind  of 
input.  An  IPA  memory  with  Fresnel  hologram  can  be  used  to  Implement  pattern 
association  with  a  large  number  of  pixels.  Experimentally,  we  have  demon¬ 
strated  it  with  24  X  24-pixel  patterns;  however,  this  by  no  means  is  the 
limit.  Because  the  resolution  of  the  holographic  mask  is  very  high,  it  does 
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not  limit  the  pixel  number.  The  maximum  number  of  pixels  is  determined  by  the 
input  mask . 


Conclusions:  An  IPA  memory  with  Fresnel  hologram  has  been  fabricated  to 
increase  the  number  of  pixels  to  be  stored  in  a  pattern.  In  addition,  the  gray 
scale  in  the  associative  matrix  is  increased  in  the  IPA  memory  to  ensure  high 
resolution  of  storage  modes  with  small  Hamming  distance. 
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[Text]  Abstract 

A  hybrid  optical-digital  system  capable  of  solving  a  set  of  linear  equations 
using  an  iterative  method  is  presented.  The  optical  part  primarily  consists  of 
a  single  holographic  mask  used  to  carry  out  vector-matrix  multiplication.  The 
remaining  part,  including  a  CCD  detector  and  a  microcomputer,  is  for  vector 
measurement  and  addition.  This  hybrid  system  was  used  to  solve  a  set  of  four¬ 
dimensional  equations.  Compared  to  theoretical  solution,  the  error  of  the 
experimental  result  is  approximately  5  percent. 

I.  Introduction 

In  recent  years,  there  are  a  number  of  studies^"^  on  optical  computation  of 
matrix  multiplication.  However,  optical  solution  of  linear  equations  is  still 
rare.  In  1982,  Casasent  presented  a  scheme  for  solving  linear  equations  by 
using  an  optical  vector-matrix  multiplier  with  acousto-optic  modulation  and 
electronic  feedback  iteration. ^  In  1986,  Caulfield,  et  al.,  presented  a  dual¬ 
mode  optical  scheme  for  solving  linear  equations  using  an  analog  and  a  digital 
optical  computing  system.^  Nevertheless,  they  failed  to  present  any  experi¬ 
mental  results  in  their  papers. 

It  is  well  known  that  there  are  numerous  iterative  methods  for  solving  linear 
equations.  The  most  common  is  a  simple  iterative  scheme.  This  scheme  involves 
the  inner  product  of  a  vector  and  a  matrix  and  a  vector  addition  in  every 
iteration.  The  former  can  be  done  optically.  As  for  the  latter,  it  is  more 
appropriate  to  use  an  electronic  system.  Therefore,  using  a  hybrid  optical- 
digital  system  to  solve  linear  equations  iteratively  is  a  very  logical  choice. 
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II .  Principle 

Let  a  set  of  linear  equations  be 


AX  -  B  (1) 

where  A  and  B  are  a  two-dimensional  matrix  and  a  one-dimensional  vector, 
respectively,  and  X  is  an  unknown  vector.  After  a  primary  transform  of  A,  the 
iterative  equation  for  this  set  of  equations  is 

X(k+i)  «  +  D  (2) 

where  matrix  T  is  related  to  A,  and  vector  D  is  a  function  of  B.  The 
superscripts  k  and  k+1  represent  the  number  of  iteration.  Equation  (2) 
converges  only  when  one  of  the  following  conditions  is  met: 

max  V  |t(i,j)  I  <  1, 

^  .  (3) 

max  |t(i,j)  \  <  1.  T  <  1, 

10<n  ifjii  ■'  j 

where  t(i,j)  is  the  element  of  matrix  T  in  the  ith  row  and  jth  column.  If  none 
of  these  three  conditions  can  be  satisfied,  it  is  necessary  to  transform  (2) 
further  to  alter  matrix  T  and  vector  B  until  t(i,j)  meets  one  of  these 
conditions . 


The  accuracy  of  solving  linear  equations  using  (2)  can  be  estimated  as 
follows : 

||X<k«)  _  x(k)||  <  e.  (4) 

where  e  is  an  arbitrary  small  quantity. 

The  hybrid  optical— digital  system  designed 
and  constructed  in  this  work  is  shown  in 
Figure  1.  In  this  system,  the  complex 
amplitude  distribution  of  the  modulating 
holographic  mask  H(m)  is  defined  by  the 
following  equations®"®; 

i  =  1,  2,  Wj,  m  =  1,  2,-,  N, 
j  =  1,  2, 

In  these  equations,  G  is  the  Fourier  transform  factor;  Ni,  N,  and  N2  are  the 
number  of  sampling  elements  on  the  input  plane,  mask  plane,  and  output  plane, 
respectively.  When  the  following  equation  is  satisfied  by  optical  parameters 

(AxAxi//A)  -  N-i,  (AxAxaZ/A)  -  N2■^  N  -=  N1N2,  (6) 
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equation  (5)  has  a  unique  solution.  Moreover,  based  on  unitarity,  it  is 
possible  to  obtain  the  following  solution  directly: 


H(m) 


t(l,j)G*  U,m) 


/sin  cir^m/U) , 


(7) 


In  equations  (6)  and  (7)  ,  Axi,  Ax,  and  Ax2  are  the  sampling  periods  of  the 
input,  holographic  mask,  and  output  phase,  respectively,  /  is  the  focal  length 
of  the  Fourier  lens,  A  is  the  wavelength,  and  is  the  ratio  of  input  plane 
sampling  element  aperture  to  sampling  period.  Therefore,  with  a  given  matrix 
(from  iteration  equation  (2)),  the  complex  amplitude  distribution  of  the 
holographic  mask  H(m)  can  be  detemined  in  equation  (7)  .  The  design  parameters 
of  the  optical  system  can  be  selected  over  a  reasonable  range  based  on 
equation  (6). 


III.  Experiment  and  Results 


The  following  is  an  arbitrary  four-element  set  of  linear  equations: 


38xi  +  3Xj  +  5X3  +  lAx^  =  70,  100X3  +  91X2  +  26X3  +  48x^  =  300,1 
9Xi  +  6X3  +  31X3  +  7X4  =  115,  X3  +  17X2  +  30X3  +  94X4  =  200,  J 

This  set  of  equations  (8)  is  transformed  into  an  iterative  format  (9): 


^X3' 

ik*l) 

'-0.24  0.06  0.10  0.28' 

fx  ^ 

(k) 

a.4> 

0.24  -0.15  0.16  0.20 

X2 

1.6 

0.18  0.12-0.38  0.14 

+ 

2.3 

^  0.01  0.17  0.30  -0.06; 

.2.0> 

Comparing  to  equation  (2),  we  get: 


0.24 

0.06 

0.10 

0.28' 

0.24  “0.15 

0.16 

0.20 

0.18 

0.12 

-0.38 

0.14 

0.01 

0.17 

0.30 

-0 . 06; 

(1.4,  1.6, 

2.3, 

2.0)T. 

(9) 


(10) 

(11) 


(T  represents  transformation.)  After  verification,  matrix  T  meets  the  conver¬ 
gence  equation  (3).  Equation  (10)  shows  that  Nj  -  N2  -  4  and  N  =  16.  Substitute 
matrix  T  into  (7)  and  let  n  “  0.4,  then  we  get  the  complex  amplitude  distribu¬ 
tion  H(m)  of  the  holographic  mask.  Its  numerical  values  are  shown  in  Table  1. 

The  design  parameters  of  the  optical  system  are  determined  based  on  equation 
(6)  and  the  existing  experimental  conditions.  They  are  AX3  =  0.2052  mm,  AX2 
0.8206  mm.  Ax  -  0.0800  mm,  /  -  415  mm,  and  A  -  6.328  x  10“*mm. 

Figure  1  shows  the  experimental  scheme  of  the  hybrid  optical-digital  system. 
The  holographic  mask  on  the  confocal  plane  was  developed  using  computer 
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generated  holograms  based  on  the  data  shown  In  Table  1.  The  amplitude  and 
phase  distribution  extend  along  the  x  axis.  Each  amplitude  phase  value  is 
represented  by  a  rectangular  grating  extending  along  the  y  direction.  If  the 
grating  period  is  p,  starting  position  is  y  and  the  light  transmitting  aper¬ 
ture  is  Ay,  then  the  first-order  diffraction  phase  of  the  grating  is  ^  -  (y/p) 
and  amplitude  A  -  [sin(7rAy/p)/(?r/p)  ]  •  li^  this  experiment,  the  grating  period 
is  0.1  mm.  The  positioning  accuracy  of  the  grating  pattern  generator  (GCA3600) 
is  0.1  fim.  Therefore,  the  phase  and  amplitude  of  the  holographic  mask  can  be 
kept  at  one  ten-thousandth.^  The  magnitude  of  the  input  vector  is  represented 
by  the  sampling  element  aperture  size  on  the  input  plane.  The  output  vector  is 
measured  by  the  CCD  detector  array.  The  result  undergoes  an  A/D  conversion  and 
then  is  input  into  a  microcomputer.  Vector  addition  is  done  in  the 
microcomputer . 

Table  1.  Amplitude-Phase  Values  of  the  Holographic  Mask  H(m) 


mm 

a 

1 

3 

4 

5 

6 

7 

8 

Amplitude 

6.211 

i.eee 

0.203 

0.193 

0.704 

0.543 

0.657 

0.622 

Phase 

1.222 

-3.057 

-1.923 

0.822 

-0.193 

-0.300 

2.033 

-0.088 

m 

9 

10 

11 

12 

13 

14 

15 

16 

Amplitude 

0.622 

0.657 

0.542 

0.704 

0.192 

0.203 

1.000 

0.311 

Phase 

O.OSS 

-2.033 

0.300 

0.193 

-0.822 

1.923 

3.057 

-1.222 

The  position  of  the  holographic  mask  is  adjusted  in  the  following  manner.  On 
the  z  axis,  the  confocal  plane  is  found  by  using  a  ground  glass  to  find  the 
maximum  spot.  A  center  line  is  set  along  the  x-axis.  It  is  lined  up  with  the 
center  line  of  sin  c  diffraction  generated  on  the  holographic  mask  plane  by 
the  input  aperture.  Alignment  along  the  y-axis  is  not  important  because  the 
grating  is  periodic. 

Let  us  assume  that  (111  1)^  is  the  initial  solution  of  the  set  of 
equations.  An  input  plate  representing  is  placed  on  plane  Xj  shown  in 
Figure  1  and  the  result  obtained  by  the  optical-digital  system  is  X^^^.  X^^^  is 
then  used  as  the  input  vector  for  the  hybrid  system  in  the  second  iteration. 
After  repeated  iterations,  the  solution  of  the  set  of  linear  equations  can  be 
obtained.  Table  2  Shows  the  experimental  data  obtained  by  the  hybrid  system 
and  the  theoretical  values  derived  from  computer-simulated  iterations.  They 
are  in  very  good  agreement. 

The  real  solution  solved  by  a  computer  is  X**  ■»  (1.01,  0.79,  3.04,  1.01)^.  After 
six  iterations,  the  solution  obtained  by  the  hybrid  optical-digital  system  has 
an  error  of 


—  y 

x'in)  -x(.n) 

2 

x’-(n) 

5.2%. 


(12) 
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Table  2.  The  Iterative  Solutions  of  the  Four-Dimensional  Equations  Produced  by 
the  Hybrid  Optical-Digital  System  and  Theoretical  Calculation 


iterative  time 

solutions  (1) 

solutions  (2) 

solutions  (3) 

solutions  (4) 

E 

■Si 

E 

WBM 

E 

T 

E 

T 

exp. 

mM 

exp. 

B3il 

exp. 

theo. 

oxp. 

tbeo 

1 

1.20 

1.20 

l.OS 

1.15  ; 

2.24 

2.24 

1.49 

1.S8 

2  1 

1.04 

0.95 

0.81 

1  0.81 

2.44 

2.57 

1.20 

1.21 

8 

0.99 

1  0.8S 

2.78 

2.84 

1.25 

1.15 

4 

0.92 

0.97 

0.77 

0.81 

2.95 

2.94 

1.11  1 

1  1.06 

5 

1.03 

0.99 

0.77 

0.80 

2.95 

2.99 

1.14  1 

1  1.03 

6 

1.05 

1.00 

0.82 

0.80 

S.02 

3.02 

1.10  1 

1.02 

Due  to  the  parallel  nature  of  the  optical  system,  linear  equations  can  be 
quickly  solved  using  a  hybrid  optical-digital  system.  However,  since  it 
involves  an  analog  operation,  its  accuracy  is  somewhat  limited.  Nevertheless, 
it  is  useful  as  the  front  end  of  the  scheme  proposed  by  Caulfield  for  using  a 
dual-mode  optical  computer  to  solve  linear  equations,  because  it  only  requires 
a  rough  solution  at  high  speed. 

The  authors  wish  to  thank  Zheng  Shihai  [6774  1597  3189]  and  Li  Dehua  [2621 
1795  5478]  for  their  assistance. 
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[Text]  Abstract 

For  the  first  time,  multiple— valued  optical  logic  has  been  demonstrated  using 
a  bistable  laser  diode  (BILD)  or  bistable  light  emitting  diode  (BILED) .  Basic 
ternary  logic  functions  such  as  COMPLEMENT,  MINIMUM,  MAXIMUM,  CYCLE,  and 
LITERALS  are  experimentally  demonstrated.  A  triplexer  circuit  design  is 
presented. 

Hybrid  optoelectronic  logic  circuits  have  attracted  a  great  deal  of  attention 
lately .  A  BILED— based  hybrid  optical  multiple-valued  logic  circuit  is 
presented  in  this  paper.  Experimentally,  five  basic  ternary  logic  functions  in 
POST  algebra,  i.e.,  COMPLEMENT,  MINIMUM,  MAXIMUM,  CYCLE,  and  LITERALS,  have 
been  implemented.  Furthermore,  a  ternary  T-gate,  or  triplexer,  design  is 
shown. 

Compared  to  binary  logic,  a  multiple— valued  logic  variable  carries  more 
information.  Therefore,  logic  operations  carried  out  by  a  multiple-valued 
logic  gate  require  less  logic  gates  and  connections.*  Just  as  Boolean  algebra 
is  the  theoretical  basis  for  binary  logic,  there  is  a  theoretical  framework 
for  multiple— valued  logic.  POST  algebra  is  a  common  one.  Considering  ternary 
logic  where  R  -  3,  its  logic  variables  are  xe{0,l,2}.  The  basic  logic 
functions  in  ternary  POST  algebra  include: 

COMPLEMENT  x  =  2  -  x, 

MINIMUM  Min  (x,y)  -=  x,  if  x  <  y. 

MAXIMUM  Max  (x,y)  -  x,  if  x  >  y. 
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CYCLE 

LITERALS 


Suc(x)  -  X  +  1,  (mod  3) 

ay.nb-f  2  If  aixib 
10  elsewhere 


The  multiple-choice  switch 
or  T— gate  operation  is: 


T(x.y.z;l0 


■I 


if  jr  =  0 
if  JC  =  1 
if  K  =  2 


Pi 


where  COMPLEMENT,  MINIMUM,  and  MAXIMUM  are  essentially  equivalent  to  NOT,  AND, 
and  OR  in  binary  logic.  CYCLE  can  be  used  to  implement  the  transformation  of 
a  logic  variable.  LITERALS  represents  the  state  of  a  logic  variable.  If  a  -  b 
-  K,  then  LITERALS  provides  the  state  of  x  at  0,  1,  2  as  Jo(x),  Ji(x) ,  and 
J2(x) .  The  T-gate  can  be  used  to  select  the  output  of  a  variable  K  at  different 
values . 

There  are  two  factors  to  be  considered 
in  implementing  multiple-valued  logic. 

(1)  Are  there  multiple  input  elements 
that  can  provide  multiple  output? 

(2)  Can  elements  be  conveniently 
connected?  It  is  simple  to  connect  ele¬ 
ments  in  an  optoelectronic  circuit.  We 
used  a  tristable  circuit  based  on  two 
BILEDs,  as  shown  in  Figure  1(a),  to 
create  a  unit  element  for  the  optoelec¬ 
tronic  ternary  logic.  In  the  figure,  the 
LED  is  a  red  GaAsP  light  emitting  diode, 

PD  is  a  3DU32  silicon  phototriode,  and  T 
is  a  transistor;  fj  and  f2  are  the  two 
optical  feedback  coefficients  for  the 
two  bistable  circuits.  The  output  ampli¬ 
tude  is  controlled  by  bias  V  and  resis¬ 
tances  Rg  and  Rj  (i  "■  1,2).  The  thresh¬ 
olds  to  turn  on  BILEDi  and  BILED2  are 

controlled  by  resistors  R^i  and  Rt2.  If  Rti  >  Rt2i  then  Poni  <  Pon2»  where 
Poni  and  Pon2  are  the  thresholds  to  turn  on  BILEDi  *nd  BILED2,  respectively.  Due 
to  the  branching  of  the  circuit,  a  tristable  output  is  achieved  with  two 
saturated  currents  and  the  switching  of  BILEDi  and  BILED2  at  different 
thresholds,  as  shown  in  Figure  1(b)  [photograph  not  reproduced]. 

If  the  three  steady  states  are  0,  1,  and  2,  and  the  switching  thresholds 
satisfy  the  following: 


Figure  1.  (a)  Schematic  Diagram  of  a 
Tristable  Circuit  or  a  Photo-Current 
Mirror;  (b)  Output-Input  Character¬ 
istic  Loop 


0  Pool  ■*'  f 
1  Pon2  2 


(1) 


then  the  circuit  can  have  a  ternary  output  with  respect  to  a  ternary  input  in 
a  reproducible  manner.  This  is  analogous  to  the  current  mirror  in  the  I^L 
multiple-valued  logic  design  in  electronics.  The  tristable  circuit  shown  in 
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Figure  1(a)  is  called  a  light-current  mirror,  a  device  which  plays  a  key  role 
in  almost  all  ternary  logic  circuits.  The  following  is  a  description  of 
optical  ternary  logic  gates  and  some  experimental  results. 


1.  COMPLEMENT  Gate 

Figure  2(a)  shows  the  circuit  of  a  ter¬ 
nary  COMPLEMENT  gate,  which  consists  of 
a  light-current  mirror  and  a  light 
emitting  diode  LED3.  Because  the  output 
of  LED3  is  opposite  to  that  of  the 
light-current  mirror,  its  output  Pg 
appears  in  a  downward  tristable  form 
with  a  triangular  wave  input  Pj .  When 
X  -  0,  the  light-current  mirror  cannot 
be  turned  on  and  all  the  current  flows 
through  LED3;  consequently,  Pg  -  2.  When 
X  -  1,  only  a  part  of  the  current  flows 
through  the  light-current  mirror,  and 
output  Pg  -  1.  When  x  -  2,  all  the 
current  flows  through  the  light-current 
mirror  and  Pg  »  0.  The  input  and  output 
waveforms  are  shown  in  Figure  2(b) 


Figure  2.  (a)  Schematic  Diagram  of  a 
Ternary  COMPLEMENT  Gate;  (b)  Input 
and  Output  Waveforms 


[photograph  not  reproduced] . 


When  the  amplitudes  are  1  and  2,  the  input  and  output  current  are  expressed  by 

the  LED  current,  which  is  9.5  mA  and  19.0  mA,  respectively.  When  Pg  -  0, 

because  of  the  added  clamping  diode  D,  the  output  of  LED3  is  essentially  held 
to  0.  At  this  time,  the  current  passing  through  LED3  is  approximately  0.8  mA. 

2.  MINIMUM  Gate 

A  ternary  MINIMUM  gate  consists  of  two  light- 
current  mirrors,  as  shown  in  Figure  3(a), 

where  BILEDi  and  BILED3,  and  BILED2  and  BILED4, 
form  two  separate  light-current  mirrors. 

Furthermore,  BILEDi  and  BILEDg,  as  well  as 
BILED3  and  BILED4 ,  are  also  connected  in 
series,  respectively.  The  switching  thresholds 
of  the  bistable  circuits  are  regulated  by 

resistors  Rt,-  (i  -  1,  2,  3,  4): 

0  <  Poni  ”  Pon2  1 
1  <  Pon3  “  ^004  ^  2  (2) 

At  this  point,  BILEDi  BILED2,  and  BILED3 

and  BILED4  are  two  separate  elements  to  be 
controlled  by  x  and  y,  respectively.  When  the 
input  signals  are  1  and  2,  they  are  switched  on  simultaneously.  If  x  and  y  are 

1,  then  BILEDi  BILED2  are  switched  on  at  the  same  time  and  the  output  is 
Pg  -  1.  If  X  and  y  are  2,  then  BILED3  and  BILED4  are  also  switched  on  and  Pg  - 

2.  Thus,  it  accomplishes  the  MINIMUM  operation  for  x  and  y.  Figure  3(b) 
[photo-graph  not  reproduced]  shows  the  experimental  results. 


(«) 

Figure  3.  (a)  Block  Diagram  of 
a  Ternary  MINIMUM  Gate; 

(b)  Experimental  Results 
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Improved  Ternary  MINIMUM  Gate 


(a) 

Figure  5.  (a)  Block  Diagram  of  a 

Ternary  MAXIMUM  Gate; 

(b)  Experimental  Results 


In  reality,  a  ternary  MINIMUM  gate  can  be  accomplished  with  an  improved  light- 
current  mirror  to  further  simplify  its  structure,  as  shown  in  Figure  4.  Two 
phototriodes  PD3  and  PD*  are  added  in  series  to  the  input  end  of  two  BILEDs. 
Since  the  photocurrent  of  the  phototriode  saturates  at  the  minimum  input  light 
intensity,  a  single  light-current  mirror  can  act  as  a  ternary  MINIMUM  gate  as 
long  as  condition  (1)  is  met. 

3.  MAXIMUM  Gate 

Figure  5(a)  shows  a  schematic  of  a  ternary  MAXIMUM  gate  comprised  of  two 
parallel-connected  light-current  mirrors.  The  light-current  mirrors  are 
similar  in  construction  to  those  described  above.  In  addition,  LED^  and  LED2, 
and  LED3  and  LED*,  are  connected  in  parallel.  Therefore,  BILEDi  and  BILEDj,  and 
BILED3  and  BILED*  form  two  separate  elements  to  complete  the  MAXIMUM  operation 
with  amplitude  greater  than  or  equal  to  1  and  that  for  amplitude  equal  to  2, 
respectively.  Once  BILEDj^  is  switched  on,  BILED2  is  also  switched  on  due  to 
light  feedback.  BILED3  and  BILED*  behave  similarly  to  BILEDi  and  BILED2.  If  the 
two  light-current  mirrors  satisfy  condition  (2) ,  then  only  BILEDi  BILED2 
are  on  when  x  —  1  or  y  *•  1,  and  Pq  “  1.  When  x  or  y  —  2,  all  the  BILEDs  are  on, 
and  Pq  -  2.  Figure  5(b)  [photograph  not  reproduced]  shows  the  experimental 
results . 

4.  CYCLE  Gate 

Figure  6(a)  is  a  schematic  of  a  one-step-right  CYCLE  gate.  A  light-current 
mirror,  comprised  of  BILEDi  and  BILED2,  is  connected  to  BILED3  in  parallel.  A 
ternary  signal  and  a  constant  light  with  amplitude  1  are  simultaneously  input 
into  the  system.  The  switching  threshold  of  BILED3  meets  the  following 
condition: 

2  <  Pon3  <  3  (3) 

Therefore,  after  the  photocurrent  jumps  up  to  2  twice  in  a  row,  BILED3  is 
switched  on  to  allow  the  current  to  flow  from  the  light-current  mirror  through 


42 


BILED3  (if  the  resistance  of  BILED3  is 
very  low  relative  to  that  of  the  light- 
current  mirror)  .  Figure  6(b)  [photograph 
not  reproduced]  shows  the  experimental 
results  of  the  one-step-right  CYCLE.  The 
input  signal  is  x  -  (210)  and  the  output 
is  Po  =  (021). 

Signal  X  can  also  be  cycled  one  step 
left,  0  2,  1  -+  0,  2  -♦  1.  Based  on 
definition,  one  step  left  is  equivalent 
to  two  steps  right.  Therefore,  one-step- 
left  may  be  accomplished  by  going 
through  the  gate  shown  in  Figure  6(a) 
twice.  Nevertheless,  it  may  also  be  done 
with  a  one-step-left  CYCLE  gate.  Figure 
7(a)  shows  a  sketch  of  the  gate.  By 

adding  BILED4  to  the  circuit  and  changing  the  constant  light  amplitude  to  2, 
the  threshold  for  switching  BILED4  on  becomes : 


Figure  6.  Schematic  Diagram  of  a 

One— Step— Right  CYCLE  Gate 


3  <  Pon4  <  4  (4) 

Hence,  after  BILED4  is  switched  on,  a  portion  of  the  current  flows  through  the 
light-current  mirror  to  result  in  Pq  -  1.  This  operation  makes  it  possible  to 
cycle  X  one  step  left. 


(a)  <*> 

Figure  7.  (a)  Block  Diagram  of  a  One-Step-Left  CYCLE  Gate; 

(b)  Circuit  Design  of  a  Ternary  Half-Adder 


It  should  also  be  pointed  out  that  if  the  constant  light  2  is  replaced  by 
another  signal  y,  the  circuit  shown  in  Figure  7  is  also  an  adder  without 
carry.  When  x+y  <  2,  the  light-current  mirror  output  is  the  direct  addition 
result.  When  x-1,  y-2,  orx-2,  y-1,  BILED3  is  turned  on  and  Pq  -  0. 
When  X  -  y  -  2,  BILED4  is  switched  on  and  Pq  -  1.  Hence,  the  operation  x  AND 
y  is  realized.  Please  note  that  when  BILED4  is  switched  on,  the  output  of 
BILED3  is  not  reduced  to  0 .  Hence ,  the  output  of  BILED3 ,  P03 ,  can  be  used  as  an 
input  for  another  BILED  where  the  threshold  is  very  low  and  the  output 
amplitude  is  1.  This  circuit  can  be  used  to  output  a  carry.  Thus,  we  have  a 
ternary  half  adder,  as  shown  in  Figure  7(b). 
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circuits  comprised  of  a  ^  ^ 

binary  NOT  gate  and  an  axcluslve  OR  gate  In  s  Multiplexer  or  T-Gate 
parallel,  as  shown  in  Figure  8.  BILEDi  and  LED4  ^ 

form  a  binary  NOT  gate,^  while  BILED2  and 

BILED3  form  an  exclusive  OR  gate.  The  output  of  LED4,  BILED2,  and  BILED3  are 
Jo(x),  Ji(x),  and  J2(x) ,  respectively.  When  x  -  0,  the  BILED  circuit  is  not 
switched  on  and  all  the  current  goes  through  LED4,  i.e.,  Jo(x)  “  2  and  Ji(x) 
=  J2(x)  “  0.  When  x  -  1,  both  BILEDi  and  BILED2  are  turned  on  and  current  flows 
through  BILEDi  and  BILED2.  Since  there  is  a  clamping  diode  in  the  LED4  branch, 
the  output  of  LED4  drops  to  0.  In  this  case,  Ji(x)  —  2  and  Jo(x)  —  J2(^)  ** 
When  X  -=  2,  the  BILED  is  also  switched  on.  Because  the  resistance  of  BILED3  is 
far  less  than  that  of  BILED2,  and  there  is  a  clamping  diode  in  the  BILED2 
circuit  as  well,  Ji(x)  falls  to  0,  i.e.,  J2(x)  -  2,  Jo(x)  -  Ji(x)  -  0. 
Experimentally,  the  output  of  BILED3  is  far  greater  than  2.  It  is  possible  to 
normalize  it  to  2  by  using  a  shunt  resistor  in  parallel  with  LED3. 


6.  T-Gate 


Let  us  re-express  the  T-gate  as  follows: 

T(x,y,z;  K)  -  Min(Jo(K),  x)  +  Min(Ji(K) ,  y)  +  Min(J2(K) ,  z)  (5) 

Therefore,  a  ternary  T— gate  consists  of  a  LITERALS  gate  and  three  MINIMUM 
gates.  Figure  9  shows  the  block  diagram  of  a  T— gate.  Through  the  LITERALS 
gate,  a  control  variable  K  results  in  three  functions  Jo(K),  Ji(K) ,  and  J2(K) . 
These  three  functions  and  the  three  signals  x,  y,  and  z  pass  through  three 
MINIMUM  gates,  respectively.  Different  variable  output  is  obtained  based  on  K. 
The  simplified  MINIMUM  gate  circuit  shown  in  Figure  4  may  be  used  as  the 
ternary  MINIMUM  gate . 
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[Text]  Abstract 

A  programmable  real-time  optoelectronic  hybrid  threshold  logic  processor  is 
presented.  Its  applications  in  optical  computing  and  optical  information 
processing  are  discussed. 

I.  Introduction 

In  optical  computing,  the  most  intriguing  subject  is  the  optical  neural  net¬ 
work,  which  mimics  the  neural  logic  processes  of  the  human  brain.  ^  Because  the 
human  brain  can  process  logic  in  a  highly  parallel,  fuzzy,  multiple— valued, 
weighted,  threshold-taking  and  feedback-controlled  manner,  it  is  necessary  to 
investigate  an  optical  logic  stronger  than  binary  vertex  primitive  logic, 
such  as  fuzzy  logic,  multiple-valued  logic,  and  binary  and  multiple-valued 
threshold  logic.  In  recent  years,  optical  multiple— valued  logic  operations 
have  been  implemented  in  several  optical  systems.^”*  Although  Arrathoon® 
suggested  using  waveguides  to  implement  optical  threshold  logic,  very  few 
simple  optical  threshold  logic  schemes  are  available.  For  the  first  time,  an 
T-Rn  (light  emitting  diode)  and  PD  (pho to triode) -based  optoelectronic  hybrid 
threshold  logic  processor  is  presented  in  this  paper.  The  processor  can  carry 
out  real-time  programmable  control  by  optical  signal  to  obtain  two  weighted 
and  threshold  values,  i.e.,  both  positive  and  negative  weighted  and  threshold 
values.  Therefore,  the  logic  processor  itself  is  a  complete  logic  set.  This 
logic  processor  is  simple  in  structure,  reliable,  and  more  importantly  com¬ 
patible  with  VLSI  (very  large-scale  integration)  technology.  Optoelectronic- 
microelectronic  integration  can  be  directly  implemented  on  a  GaAs  or  GaAs/Si 
wafer.  Hence,  it  has  a  bright  future  in  optical  computing  and  digital  optical 
information  processing. 
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II.  Optoelectronic  Hybrid  Threshold  Logic 
Processor 

Figure  1  shows  the  schematic  diagram  of  a 
threshold  logic  gate.  It  has  n  binary  input 

variables  Xf(i-1,2, . ,n).  Each  variable 

is  accompanied  by  a  weight  factor 

Wi(i-1,2, . ,n) .  There  is  a  relation 

between  output  y  and  expected  threshold  T®: 


1, 

0, 

WjXi<T. 

Only  when  the  wei^ted  input  variable  is 
greater  than  or  equal  to  the  threshold  is 


Figure  1.  Schematic  Diagram  of 
Threshold  Gate 


the  output  y  eqtial  to  1.  Otherwise,  it  is  0.  If  a  switching  function  /(x)  is 
a  threshold  logic  function,  i.e.,  it  can  be  produced  by  a  threshold  logic 
gate,  then  there  is  a  weight-threshold  set  {wi,W2,  •  •  •  ,Wn;T)  that  satisfies 
equation  (1).  Mathematically,  if  {w,-,T)  covers  all  real  numbers,  then  a  single 
threshold  logic  gate  forms  a  complete  logic  set.  This  means  that  all  logic 
functions  can  be  implemented  by  a  finite  number  of  threshold  logic  gates. 
Furthermore,  a  threshold  logic  gate  has  more  input  capacity  than  a  binary 
basic  logic  gate.  Hence  the  number  of  logic  elements  and  transmission  lines 
can  be  significantly  reduced. 


Figure  2.  (a)  Circuit  of  Optoelectronics  Threshold  Logic 

Processor;  (b)  Part  of  the  Circuit  of  Figure  2 
(a)  Which  Only  Consists  of  Input  X{ 


Figure  2(a)  shows  the  circuit  of  an  optoelectronic  hybrid  threshold  logic 
processor.  In  reality,  it  is  a  variation  of  a  bistable  LED  circuit.*  Experi¬ 
mentally,  optical  feedback’  was  removed  to  increase  the  response  speed  of  the 
logic  gate.  The  LED  output  is  the  logic  gate  output  y.  Each  weighted  input  end 
contains  four  silicon  PDs  in  series.  Two  serve  as  the  input  of  the  input 
variable  and  the  other  two  are  the  input  of  the  weight  signal.  Each  weij^t  Wf 
is  divided  into  two  parts:  Wj"*"  and  w,-".  When  Xf  is  positively  weighted. 
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Wj  -  and  w,-~  -  0.  If  x,-  is  negatively  weighted,  then  w,-  -  -w,-"  and  w,-"^  =  0. 
Threshold  T  is  also  controlled  by  two  PDs  connected  in  series.  Similarly,  it 
has  a  positive  and  negative  value.  Different  weight  and  threshold  values  can 
be  obtained  by  varying  the  intensities  of  the  weight  and  threshold  signal. 
Consequently,  this  threshold  logic  processor  becomes  real-time  programmable. 

In  order  to  illustrate  how  real-time  weight  and  threshold  can  be  implemented 
by  the  processor,  let  us  use  an  example  where  there  is  only  one  input  variable 
Xi ,  as  shown  in  Figure  2(b).  The  input  Xf  lHuuninates  PD,-'*’  and  PDj”  and  the 
weight  control  signals  w,-'''  and  w,-”  shine  on  and  ,  respectively.  The 

current  intensity  Iji  can  be  expressed  as: 

Asstiming  all  PDs  are  identical  and  operate  in  the  linear  amplification  area, 
then  the  output  photocurrent  Ij  is  proportional  to  the  input  light  intensity 
P{,  i.e. ,  Iji  -  KP,*,  where  K  is  an  optoelectronic  conversion  constant.  Becausepo] 
and  pD*^  (PD,“  and  pD-^)  are  connected  in  series,  (jj^)  is  proportional  to 
either  the  x^  or  Wj'*’  (wj") ,  whichever  has  the  lower  optical  intensity;  i.e.: 

=  K min (P(Xi) ,  Piw*))  (3) 

where  min(x)  is  the  minimxam  threshold  function.  P(xj)  and  p(Wi)  are  the 
intensities  of  x,-  and  w* ,  respectively  .  The  way  weight  and  threshold  values  are 
changed  is  different  from  multiplying  the  amplitude  of  the  logic  variable  by 
a  factor  w,-.  From  equation  (3)  we  can  see  that  p{w±)  Is  far  less  than  the 
intensity  at  xj  —  1.  Hence,  is  a  function  of  p{wi)  only  when  x,-  -  1.  We 
have : 

_*  _  jKP(Wi) ,  for  Xi  *  1  (4) 

0,  forx.^O  ^  ^ 

Since  P(Xi)  |;,^.i>P(iir|)  ,  let  PU^)  =  gP(wl)  |„..i ;  this  means  that  the 

intensity  amplitude  at  x,-  -  1  is  evenly  divided  into  q  parts.  If  the  weight  of 
Xj  is  ±r(l  ^  r  ^  q),  then  p(wi)  =  rP(Wi)  |„..i .  For  example,  at  xj  -  1,  P(x,-)  - 
50  pW  and  at  wl  -  1.  P(Wi)  -  5  pW;  then  q  -  10.  If  the  weight  of  x,-  is  5,  then 
the  intensity  for  w,"^  is  25  pM.  Of  course,  the  value  of  q  may  be  chosen 
arbitrarily.  However,  P(Wi)  |„..i  must  be  greater  than  the  saturated  threshold 
of  the  PD.  Usually,  this  threshold  is  very  low,  approximately  0.5  pW.  Hence, 
the  weight  may  vary  over  a  wide  range. 

Let  us  assume  that  there  are  m  positive-weight  input  signals  and  n-m  negative- 
weight  input  signals  in  the  circuit  of  Figure  2(a).  Then,  the  current  after 
taking  the  threshold  is 

KminlPiXi),  P(wl)}  Kmin{P{Xj),  P(wj)]  ±  AP(D  (5) 
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Here,  T  -  when  ±KP(T)  is  negative,  and  T  -  -T"  when  it  is  positive. 
Considering  the  effect  of  minimum  saturated  threshold  intensity  described 
above,  the  threshold  intensity  P(T)  is  adjusted  so  that  P(T^)  -  P(T)  -  0.5  pW 
P(T~)  “  P(T)  +0.5  pW.  This  is  equivalent  to  adding  a  constant  maintenance 
signal  H  where  P(H)  «  0.5  pW  on  top  of  PDt“  to  overcome  the  saturated  threshold 
light  intensity.  In  this  case,  the  threshold  logic  output  is  y  -  1  when 
Id  >  0.  Different  logic  functions  can  be  obtained  by  choosing  different  weight 
and  threshold  values. 


Figure  3.  Experimental  Results  for  Input  Variables  Nvimber  n  -  3. 

The  upper  three  waveforms  Indicate  Xi,  X2,  and  X3,  respec¬ 
tively;  the  bottom  waveform  shows  the  output  Y. 

(a)  AND;  (b)  OR;  (c)  MAJORITY;  (d)  /(x)  -  Xj  +  X2X3; 

(e)  /(x)  -  X1X3  +  X2(xi+X3);  (f)  /(x)  -  x^  +  X2  +  Xg 

Figures  3  (a)  -  (f)  are  the  experimental  results  with  n  -  3.  The  three  signals 
were  supplied  by  three  LEDs  driven  by  three  s3nichronized  square-wave  voltage 
signals  produced  directly  by  the  signal  generator.  The  peak  voltage  is  10  V 
and  the  current-limiting  resistor  has  a  value  of  220  0.  The  weight  signalwj 
and  threshold  T  are  also  supplied  by  LEDs.  By  using  switches  and  variable 
resistors,  it  is  possible  to  change  their  signs  and  amplitudes.  The  voltage 
for  the  driver  and  the  main  circuit  is  Vcc  -  4.5  V.  The  series  resistance 
values  corresponding  to  w,-  -  1  and  2  are  1.2  kO  and  600  0,  respectively.  Figure 
3(a)  and  (b)  shows  the  AND  and  OR  logic  at  T  -  1  and  3  when  w^  -  1.  Figure  3(c) 
shows  the  MAJORITY  logic  when  w^  —  1  and  T  —  2.  When  more  than  two  of  Xi,  Xg, 
and  X3  are  1,  the  output  is  1.  Similarly,  when  wj  “  —1  and  T  —2,  we  have  the 
MINORITY  logic.  Obviously,  the  only  way  to  get  a  more  abundant  logic  function 
is  to  have  different  input  weight  values .  For  instance ,  the  weight- threshold 
set  {2, -1,1;!}  can  result  in  /(x)  -  Xj  +  jc^Xg.  With  a  set  such  as  {1,-1, 1;1}, 
we  have  /(x)  -  XjXg  +  X2(xi+X3).  With  {-1,-1, 1;-1)  we  have  /(x)  -  Xj  +  Xg  +  Xg. 
These  results  are  shown  in  Figures  3(d)-(f),  and  represent  only  a  small 
portion  of  the  logic  functions  attainable  from  a  single  threshold  logic  gate. 
It  should  be  pointed  out  that  only  six  out  of  eight  of  the  possible 
combinations  of  the  three  input  signals  are  shown  in  Figures  3(a)-(f).  Two 
combinations,  (0,1,0)  and  (1,0,1),  are  not  included.  However,  this  input 
signal  arrangement  can  provide  the  greatest  number  of  combinations  possible. 
Therefore,  Figures  3(a)-(f)  can  be  considered  to  reflect  the  true  values  of 
logic  operation. 
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A  threshold  logic  gate  can  have  multiple 
binary  output  values  such  as  y^,  72,  •  •  • , 
y,„.  Under  identical  weighting,  each  out¬ 
put  end  provides  a  different  threshold. 

This  requires  the  reproduction  of  the 
weight  into  m  parts  which  should  not 
interfere  with  each  other  when  the 
threshold  is  applied.  Figure  4  shows  a 
logic  processor  with  double  output.  The 
processor  has  two  amplifier  FDs  with 
opposite  output  characteristics;  one  is 
an  NPN  triode  Ti  and  the  other  is  a  PNP 
triode  T2.  LEDi  and  LED2  are  the  output 
ends.  -Vcc  is  the  negative  bias.  The  out¬ 
put  of  LED2  may  be  1  only  with  a  nega¬ 
tive  weight  input  and  greater  than 
negative  threshold  value.  In  principle, 

n 

when  ^  WjXi  =  T,  both  outputs  yi  and  72 

should  be  1.  In  this  case,  if  the  weight-threshold  set  corresponding  to  yi  is 
{Wi,W2, • • • ,Wn;T} ,  then  the  weight-threshold  set  corresponding  to  72  is 
{-Wi,-W2,  •  •  •  ,-Wn;-T) .  However,  because  of  the  cutoff  property  of  Tj  and  T2,  also 
because  it  cannot  make  correction  for  a  threshold  simultaneously,  it  is  only 
possible  to  add  H  on  the  input  ends  T“  and  T*"  for  correction.  Therefore,  in  a 
simple  manner  yi  and  are  opposite  in  nature,  i.e. ,  yi  -  72.  If  the  threshold 
bias  maintenance  light  H  is  removed, it  requires  another  state  yi  ”  72  “  0.  In 
the  following  discussion,  the  applications  of  such  logic  devices  will  be 
analyzed. 

III.  Discussion 

Threshold  logic  is  not  only  an  effective  means  to  implement  the  logic 
functions  in  multi-variable  Boolean  algebra,  but  also  can  be  widely  used  in 
digital  graphic  processing,  mode  recognition,  classification  of  fuzzy 
functions,  and  neural  networks.  Existing  optical  neural  network  models,  such 
as  the  Hopfield  model, ^  are  essentially  based  on  threshold  logic.  Simple 
processes  such  as  memory  connection  and  learning  are  completed  by  a  weighted 
feedback  input  and  threshold  output  between  neurons.®  Therefore,  threshold 
logic  can  be  directly  applied  in  neural  network  research.  As  a  matter  of  fact. 
Figure  2(b)  can  be  considered  as  a  simple  optoelectronic  hybrid  neural 
element.  It  can  directly  add  and  subtract  weighted  EXCITATORY  input  and 
weighted  INHIBITORY  input.  Hence,  this  neural  element  can  perform  as  a 
threshold  gate.  In  comparison,  the  threshold  gate  in  existing  neural  networks 
requires  very  complicated  electronics  and  a  computer.  In  addition,  this  neural 
element  can  have  bipolar  binary  outputs  of  +1  and  -1  (see  Figure  4) .  VLSI 
technology  and  free-space  optical  connection  technology  will  make  this  type  of 
neural  network  even  more  attractive. 


Figure  4.  Schematic  Diagram  of 
Threshold  Logic  Processor  With 
Double  Output  Y^  and  Y2,  Where 
Yi  -  /(X)  and  Ya  -  /(x) 
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Two-dimensional  parallel  logic  processing  is  also  a  major  subject  in  optical 
computing.®  Optical  threshold  logic  also  has  some  potential  in  this  area.  After 
connecting  an  OR  gate  in  series  with  the  threshold  logic  processor  for  the 
double  variable  input  shown  in  Figure  4,  16  Boolean  logic  functions  can  be 
implemented.  In  general,  can  accomplish  14  logic  functions  other  than  XOR 
[i.e.,  exclusive  OR]  and  XNOR  [i.e.,  exclusive  NOR].  If  threshold  bias  is 
removed  from  the  circuit  of  Figure  4,  I.e. ,  yj  -  yj  -  0  is  allowed  to  exist, 
then  outputs  yi  and  yg  represent  Xi-5c^  and  Xi-X2,  respectively,  when  Wj  -  1, 

W2  -  -1,  and  T  -  0.  When  an  OR  gate  level  is  added  in,  the  OR  gate's  output 
y  -  yi  +  y2  directly  yields  the  XOR  logic.  Similarly,  when  Wi  -  W2  -  1  and 
T  —  1,  y  directly  provides  XNOR  logic.  Such  a  logic  processor  can  serve  as  a 
pixel  in  a  two-dimensional  parallel  logic  process.  It  can  be  directly 
controlled  by  an  optical  signal  to  Implement  real-time  programmable  space- 
invariant  and  space-variant  logic. 
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